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Introduction 


The  principal  objective  of  this  proposal  is  to  obtain  physico-chemical  information  on  estrogen 
derivatives  and  correlate  this  knowledge  to  their  knovra  biological  functionality.  Slight  chemical 
variation  in  these  molecules  can  change  the  carcinostatic  potentials  from  agonistic  to  inhibitory. 
Relating  these  biological  reactions  to  physical  properties  such  as  point  charges  of  atoms  and  the 
electrostatic  potential  is  a  logical  first  step  in  the  intelligent  design  of  therapeutical  drugs.  We  are 
obtaining  information  about  the  electronic  properties  of  estrogen  derivatives  from  experimental 
determination  of  their  electron  density  using  high  quality  single  crystal  X-ray  crystallography  at 
low  temperatures. 


Body 

Task  1.  Preliminary  Studies  on  a  series  of  crystals  of  estrogen  derivatives 

•  Feasibility  study  of  already  available  crystals  of  estrogen  derivatives 

The  2,3,17-estriol  crystals  appeared  to  be  air  and  light  sensitive.  We  are  planning  to  work 
on  crystallization  of  this  compound  under  nitrogen  atmosphere. 

•  Development  of  crystallization  methods  for  the  derivatives  not  vet  available  as  high 
quality  single  crystals 

We  have  continued  to  develop  crystallization  methods  for  obtaining  high  quality  crystals 
of  several  different  estrogen  derivatives.  During  the  past  year  we  succeeded  to  obtain 
high  quality  crystals  of 

1 7a-estradiol  (crystallized  from  ethanol), 

-  4-Hydroxyestradiol  (crystallized  from  acetonitrile) 


•  Variable  Temperature  studies  on  each  sample  to  define  the  appropriate  temperature  for 
the  measurement 


Variable  temperature  studies  were  performed  on  the  diffractometer  on  both  above 
mentioned  crystals.  As  our  criterion  for  stability  we  chose  the  change  in  mosaicity  of  the 
respective  crystal  as  identified  by  the  broadening  of  the  profiles  of  the  diffraction 
intensities  in  comparison  to  a  room  temperature  measurement.  In  our  in-house 
experiments  using  Ag-radiation  (0.56083  A),  both  samples  proved  to  be  stable  under 
nitrogen  cooling  to  100  K. 

The  3-Hydroxy-P-estradiol«MeOH,  3-Hydroxy-P-estradiol*urea«MeOH»H20,  3- 
Hydroxy-a-estradiol*MeOH,  Hydroxy-a-estradiol«urea*MeOH*H20,  P-estradiol 
•isopropanol,  estrone»acetonitrile,  estriol»acetone  crystals  were  tested  under  helium 
temperature  (about  15K).  Contrary  to  our  previous  synchrotron  results  with  17  a- 
estradiol«l/2  H2O,  17P-estradiol«urea,  17P-estradiol«l/2  MeOH,  17P-estradiol»2/3 
MeOH*l/3  H2O  and  estrone,  we  found  no  significant  increase  in  the  crystal's  mosaicities 
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after  slow  cooling.  We  are  continuing  our  investigations  of  the  correlation  between  the 
final  temperature  of  the  crystal,  speed  of  cooling,  flux  of  the  beam,  wavelength  etc.  and 
the  mosaicity  observed  from  the  diffracted  beam  in  order  to  enhance  the  stability  of  the 
respective  estrogen  derivative. 

•  Routine  X-ray  crystal  structure  determination  on  uncharacterized  derivatives 

crystal  structures  of  all  studied  compounds  have  been  reported  earlier 
Task2.  Electron  density  studies  on  the  above  mentioned  estrogen  analogues 

•  Electron  density  studies  at  nitrogen  temperatures 

We  have  collected  complete  high  quality  data  sets  for  charge  density  analyses  on  crystals 
of 

1 7a-estradiol 

-  3,16a,17P-estriol 
17p-estradiol*l/2  H2O 

-  1 7p-estradiol»urea-  two  data  sets  were  collected  due  to  reflection  overlapping 
problems  in  the  first  data  set. 


•  Electron  density  studies  at  helium  temperatures 

He-temperature  studies  have  shown  contradictory  results  with  respect  to  crystal 
instability  (see  above)  under  extreme  cooling.  We  are  developing  cooling  techniques  that 
will  improve  the  stability  of  the  crystals  and  therefore  enable  us  to  use  the  benefit  of  very 
low  temperatures  (a  vast  increase  in  scattering  power)  for  our  electron  density  data  sets. 

•  Analysis  of  the  experimental  data,  preparation  of  manuscripts 

-  3,16a,17P-estriol;  The  electron  density  study  is  completed,  a  manuscript  for 
publication  is  in  preparation.  An  optimized  local  atomic  coordinate  system  was 
developed  (Appendix  A)  in  order  to  better  compare  the  results  from  different 
estrogen  analogs.  The  results  (Appendix  B)  were  presented  at  the  Gordon 
Research  Conference  "Electron  distribution  and  chemical  bonding",  Mount 
Holyoke  College  -  South  Hadley  ,  MA  July  8-  12, 2001. 

Estrone:  A  manuscript  for  publication  has  been  prepared  and  is  ready  for 
submission  (Appendix  C). 

-  4-Hydroxyestradiol:  Preliminary  crystallographic  studies  have  been  performed, 
the  crystals  appeared  to  be  good  quality,  the  crystal  structure  (Appendix  D), 
redetermined  at  120K  confirmed  the  results  reported  earlier  from  room 
temperature  data. 

-  17a-estradiol«l/2  H2O:  The  analysis  of  the  electron  density  data  is  currently 
under  investigation  but  incomplete. 
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17a-estradiol:  The  analysis  of  the  electron  density  data  is  currently  under 
investigation  but  incomplete. 

17P-estradiol*l/2  MeOH:  The  analysis  of  the  electron  density  data  is  currently 
under  investigation  but  incomplete. 

17P-estradiol*urea:  The  analysis  of  the  electron  density  data  is  currently  under 
investigation  but  incomplete. 

-  17p-estradiol*2/3  MeOH*  1/3  H2O:  The  publication  of  the  crystal  structure  is  in 
preparation.  Analyzing  the  electron  density  data  has  proven  to  be  difficult, 
because  the  compound  crystallizes  in  the  space  group  P2i  with  three 
crystallographically  independent  estradiol  molecules  in  the  asymmetric  imit.  The 
number  of  parameters  therefore  triples  compared  to  crystal  structures  with  only 
one  estrogen  molecule  in  the  asymmetric  unit.  We  believe  that  the  results  of 
electron  density  determinations  on  similar  derivatives  will  provide  us  with 
valuable  starting  values  to  further  model  the  electronic  properties  of  the  crystal 
under  investigation. 

17P-estradiol*l/2  H2O:  The  analysis  of  the  electron  density  data,  measured  at  90K 
with  Bruker's  CCD-6000  detector  showed  an  unsatisfactory  quality  of  the  data 
due  to  an  unknown  measurement/data  integration  problem.  Currently,  we  are 
trying  to  understand  the  source  of  the  poor  quality  of  the  data  coming  from  the 
measurements  with  CCD-6000.  Our  preliminary  recommendation  is  not  to  use 
this  detector  model  for  accurate  charge  density  studies,  at  least  in  combination 
with  a  Ag  X-ray  tube.  The  data  will  to  be  recollected  with  CCD  2K. 

Task  3.  Extension  of  the  above  mentioned  studies  to  the  series  of  E?  C-rine  analogs. 

The  compounds  have  been  ordered  and  we  are  planning  to  start  the  feasibility  tests  and 

try  to  crystallize  these  samples  shortly. 


Key  Research  Accomplishments 

•  Crystallization  of  high  quality  crystals 

•  Determination  of  electrostatic  potential  and  other  properties  of  3,16a,17p-estriol 
and  estrone 

•  Methodology  development  for  electron  density  experiments  and  data  treatment: 

a)  optimized  temperature  conditions 

b)  best  detector  type  determination 

c)  best  local  atomic  coordinate  system 


Reportable  Outcomes 


1.  American  Crystallographic  Association  Meeting,  St.  Paul,  MN,  July  22-17,  2000. 

2.  1 9th  European  Crystallographic  Meeting,  Nancy,  France,  August  25-31,  2000. 

3.  Gordon  research  conference  "Eleetron  distribution  and  chemical  bonding".  Mount 
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Holyoke  College  -  South  Hadley  ,  MA  July  8  -  12, 2001. 

4.  Bruker  Area  Detector  Users  Group  Meeting,  Madison,  WI,  May  21-22, 2001. 


Conclusion 


We  have  derived  the  electron  density,  electrostatic  potential  and  related  properties  for 
two  estrogen  derivatives:  3,16a,17P-estriol  and  estrone.  Though  it  is  too  early  to  make  any 
predictions,  these  two  compounds  show  a  significant  difference  in  the  electrostatic  potential 
distributions  and  atomic  charges.  We  have  developed  the  methodology  of  crystallization  and  X- 
ray  experiments  to  obtain  even  more  accurate  results  in  a  shorter  period  of  time.  For  four  other 
estrogen  derivatives  charge  density  quality  data  have  been  collected.  These  data  are  currently 
under  multipole  refinements. 
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Appendix  A. 


An  Optimized  Atomic  Local  Coordinate  System  for  Estrogen  Molecules 


Charge  Density  'Rules’: 

1 .  Non-hydrogen  atoms,  always  label  the  first  vector  (vl)  x.  (unless  it's  a  linear  atom,  then  vl 
should  be  z). 

2.  Always  label  the  second  vector  (v2)  y. 

3 .  Hydrogen  atoms,  always  label  the  first  vector  (v  1 )  z. 

What  to  populate  first: 

sp  20,  possibly  40  (watch  pi-density) 

sp2  33+,  20 

sp3  32-  bisect  2  pairs  of  atoms  for  dummy  atom  positions 
Remember:  always  make  the  bond  distance  of  the  atoms  equal. 

For  example,  if  you  want  to  bisect  the  angle  of  a  H  atom  and  C  atom,  push  the  H  to  the  C 
distance,  or  your  vector  will  not  be  a  true  bisection. 


Treat  lone  pairs  like  bonds  when  determining  hybridization. 
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)p  side,  or  side  containing  the  methyl  group,  is  the  p  side.  Therefor  all  hydrogens  on 
de  are  to  be  named  with  a  “b”.  All  hydrogens  on  the  bottom  side,  or  a  side,  will  be 
d  with  an  “a”. 
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Che  correct  enantiomer,  ii 
the  angle  is  positive  (+), 
you  have  the  wrong 

enantiomer.  BOTTOM 


M 


R2CH2  Hydrogens 

These  include  the  hydrogens  on  the  following  C  atoms. 

(C6,C7,C11,C12,C15,C16) 

This  example  is  of  C6. 


Point  v2  at  the 
other  hydrogen 
on  the  carbon. 


11 


H8,  H9,  and  H14 


H9  C7 
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Hydroxy  Group 

This  diagram  can  be  used  for  any  hydroxy  group  (X=  C,  N,  etc,) 


To  calculate  the  bisection, 
place  a  hydrogen  atom  at  the 
lone  pair  position,  then  push  it 
to  the  same  distance  as  the 
O-X  bond. 

Remember  to  pay  attention  to 
the  geometry.  Sometimes  the 
oxygens  can  appear  to  be 
shaped  more  like  sp3  than  sp2. 
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C 1 8  with  Hydrogens 


C13  v2 
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16 


HI  la'" . 


CIO^  / 
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This  diagram  is  for  17p-estradiol. 

If  you  have  a  different  derivative, 
C16&  Cl 7  might  differ.  IfC16& 
C17  remain  sp  the  current  scheme 
will  work  by  simply  substituting  the 
different  atoms.  If  they  are  a 
different  hybridization,  scheme 
must  change  to  support  the 
appropriate  multipoles. 


C8, 


D 

Ring 


:.ci4 
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A  Ring  Hydrogens 


19 


Coordinate  System  Atom  List 

itoms  in  parenthesis  (),  the  vector  should  point  to  a  dummy  atom  bisecting  these  2 


If  there  are  2  atoms  in  parenthesis 
atoms 


Cl 

C2 

CIO 

C2 

C3 

Cl 

C3 

C4 

C2 

C4 

C5 

C3 

C5 

CIO 

C4 

C6 

(C5-C7) 

(C5-H6a) 

C7 

(C6-C8) 

(C6-H7a) 

C8 

(C7-C9) 

(C7-C14) 

C9 

(C8-C10) 

(C8-C11) 

CIO 

Cl 

C5 

Cll 

(C9-C12) 

(C12-Hlla) 

(C11-C13) 

(C13-H12a) 

C13 

(C12-C14) 

(C14-C17) 

CM 

(C8-C13) 

(C8-C15) 

C15 

(C14-C16) 

(C14-H15a) 

C16 

(C15-C17) 

(C15-H16a) 

C17 

(C13-C16) 

(Cl  6-02) 

C18 

(H18a-H18b) 

(H18b-H18c) 

Hvdroaen  Neutron 


Distances 


-CH3  =  1.059 


-C2-CH2=  1.092 


-Car-H  =  1.083 


-C3-CH=  1.099 


-C-OH  =  0.967 


HI 

Cl 

C9 

H2 

C2 

HI 

H4 

C4 

01 

H6a 

C6 

H6b 

H6b 

C6 

H6a 

H7a 

C7 

H7b 

H7b 

C7 

H7a 

H8 

C8 

C13 

H9 

C9 

H14 

Hlla 

Cll 

Hllb 

Hllb 

Cll 

Hlla 

H12a 

C12 

H12b 

H12b 

C12 

H12a 

H14 

CM 

H9 

H15a 

C15 

H15b 

H15b 

C15 

H15a 

H16a 

C16 

H16b 

H16b 

C16 

H16a 

H17a 

C17 

02 

H18a 

C18 

H18b 

H18b 

C18 

H18c 

H18c 

C18 

H18a 

% 
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A  Charge  Density  and  Electrostatic  Potential  Study  on  Estrogen  Molecules 

16a,  17P-Estriol 
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A  charge  density  and  electrostatic  potential  study  on  estrogen  molecules: 

16a,  17(3  -  Estriol 

Abstract 

The  measurement  of  the  electron  density  distribution  of  the  principal  members  of  the 
estrogen  family,  the  derivatives  estriol,  estradiol  and  estrone  are  in  progress.  As  a  part  of  our 
study  the  electron  density  distribution  of  16a,17P-estriol  (C18H24O3)  was  determined  from  high 
resolution  single  crystal  x-ray  diffraction  measurements  at  100  K.  Estriol  crystallizes  with  two 
molecules  in  the  asymmetric  unit.  Intensity  data  were  collected  using  a  SMART  2K  CCD  area 
detector  with  Ag  Ka  radiation  to  a  resolution  corresponding  to  (sin0/X,)max=l-  29  A-'.  The 
structure  was  solved  and  a  conventional  spherical  atom  refinement  was  carried  out.  A  multipole 
pseudo  atom  refinement  was  then  performed  using  the  XD  program.  The  electron  density 
multipolar  parameters  were  used  to  calculate  the  electrostatic  potential  around  the  molecule 
removed  from  the  crystal  lattice.  Strong  electro  negative  regions  were  found  in  the  vicinity  of  the 
oxygen  atoms  in  both  the  molecules.  Net  atomic  charges  were  calculated  from  a  kappa 
refinement.  The  topological  analysis  of  the  hydrogen  bond  indicates  the  strength  of  the 
intermolecular  interactions. 

Introduction 

More  than  60%  of  all  breast  cancers  are  known  to  be  hormone  dependent,  ie.,  initiation 
and  progress  can  be  influenced  by  estrogens  and  related  compounds.  These  molecules  bind  as 
ligands  to  the  estrogen  receptor  and  hereby  initiate  a  series  of  events  resulting  in  the  activation  or 
repression  of  selective  genes.  Subtle  changes  in  the  structure  of  the  estrogen  affect  their 
chemical/biological  behavior  and  can  result  in  the  development  of  growth  inhibitors  for  tumors. 
For  example  moving  the  hydroxy  group  from  C(3)  (estradiol)  to  one  of  the  adjacent  C-atoms  can 
change  the  carcinostatic  potentials  from  agonistic  (estradiol)  to  inhibitory  (  2-and  4-hydroxy 
estratrien  -17p-ol).  While  most  changes  modify  the  affinity  of  the  ligand  for  the  estrogen 
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receptor,  this  does  not  necessarily  correlate  with  the  stimulation  ability  of  the  transcription  of 
estrogen  responsive  genes.  The  mechanism  by  which  the  ligands  regulate  the  gene  expression  is 
currently  unknown,  however  it  has  been  suggested  that  the  differences  in  electrostatic  potential 
are  responsible  for  the  variation  in  the  regulation  of  hormone  dependent  genes.  This  information 
can  be  obtained  from  studies  of  the  electronic  properties  of  estrogen  derivatives.  Experimental 
determination  of  the  electron  density  distribution  of  molecules  by  x-ray  diffraction  methods 
leads  directly  to  the  determination  of  properties  such  as  charges  of  atoms  and  the  electrostatic 
potential,  which  in  turn  aid  in  the  prediction  of  how  this  molecule  might  react  with  a  given 
environment. 

We  have  begun  a  systematic  study  of  such  compounds  by  investigating  crystals  of  the 
derivatives  of  the  principal  estrogens,  estriol,  estradiol  and  estrone.  The  second  study  in  this 
series,  on  estra-l,3,5(10)-triene-3,16a,17p-triol,  is  reported  here.  Estriol  is  the  one  of  the  active 
estrogens  found  in  the  body.  The  main  source  of  estriol  is  the  placenta,  which  produces  a  large 
amount  of  estriol  during  pregnancy.  The  commercially  available  estriol  has  been  crystallized 
from  acetonitrile.  The  geometry  and  electronic  structure  of  the  two  molecules  in  the  asymmetric 
unit  as  described  below. 


Experimental  Section 


1.  Data  collection: 

The  estriol  was  crystallized  from  acetonitrile  by  slow  evaporation  at  room  temperature. 
The  crystals  were  found  to  be  colorless  and  prismatic.  In  order  to  find  a  sample  with  good 
reflection  profiles  several  crystals  were  examined  on  the  diffractometer  for  low  and  high  angle 
scattering.  A  well  shaped  0.48  x  0.26  x  0.24  mm  size  crystal  was  selected  and  mounted  on  a 
Bruker  SMART  2K  CCD  diffractometer  using  graphite  monochromatized  Ag  Ka  (A,=0.56089A) 
radiation  (50V,  35mA).  As  a  preliminary  check  a  room  temperature  data  set  was  collected  for  a 
hemisphere.  Refinement  of  the  structure  from  this  data  gave  results  comparable  to  the  reported 
structure.  Then  the  crystal  was  cooled  to  100  K  using  an  Oxford  Cryostream  N2  open  flow 
cryostat.  When  the  crystal  was  cooled,  no  significant  unit  cell  modifieation  was  observed  and  it 
was  very  stable  over  the  period  of  the  entire  data  collection. 
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Initially  the  orientation  matrix  and  the  unit  cell  parameters  were  obtained  from  a  least- 
squares  refinement  using  88  well  centred  reflections  measured  in  three  sets  of  matrix  runs  of 
each  20  frames.  The  low  temperature  cell  thus  obtained  is  ~2  %  smaller  than  that  at  room 
temperature.  A  full  sphere  of  low  and  high  angle  reflection  intensities  were  measured  using  two 
detector  settings  at  20  =  0°  and  20=  -50°,  using  the  (o  scan  method  for  different  (p  angles  (  0, 45, 
90, 135,  180, 225, 270  and  315°).  The  CCD  detector  was  placed  at  5.1  cm  from  the  crystal  where 
the  reflection  overlap  was  minimal.  A  total  of  600  frames  were  collected  for  each  cp  angle  with 
a  scan  width  of  0.3°.  The  frame  times  for  low  (40  sec)  and  high  angles  (160  sec)  were  decided 
after  careful  observation  of  the  most  intense  reflection  thus  avoiding  detector  overflow  but 
utilizing  the  full  dynamic  range.  The  crystal  decay  was  monitored  by  re-measuring  the  initial  50 
frames  at  the  end  of  each  data  collection,  and  it  was  found  to  be  negligible.  The  data  collection 
strategy  was  checked  with  the  ASTRO  program  incorporated  into  the  SMART  software.  We 
obtained  100%  coverage  for  a  resolution  of  (sin0/X)max  =  1-1  A''  with  good  redundancy.  A  total 
of  9600  frames  were  collected  over  the  period  of  two  weeks. 

2.  Data  reduction  and  Averaging 

The  best  unit  cell  and  the  orientation  matrix  were  determined  by  thresholding  [I/a(I)  > 
20]  of  reflections  from  each  run  of  the  data  set.  This  orientation  matrix  fits  well  with  the  entire 
reciprocal  space  observed.  The  data  has  been  integrated  and  corrected  for  background,  Lorentz 
and  polarization  effects  using  SAINT  sofware's  advanced  integrating  procedures.  The  method  of 
integration  of  the  collected  frames  was  as  described  by  Kabch  (1993)  for  area  detectors.  During 
the  integration  the  orientation  matrix  and  unit  cell  were  not  allowed  to  vary.  Appropriate  three 
dimensional  box  sizes  were  used  in  low  and  high  angle  integrations  to  account  for  the  Kai  and 
Ka2  splitting  in  high  order  reflections.  No  absorption  correction  was  applied  because  of  the  low 
value  of  the  absorption  coefficient  (  0.055  mm"’).  All  the  corrected  reflection  intensities  were 
scaled  and  averaged  based  on  the  point  group  symmetry  using  the  program  SORTAV.  A  total  of 
278184  reflections  were  measured  to  sin0/A,  =  1.1  A'*.  33019  are  unique  reflections  with  an 
average  redundancy  of  7.3.  Of  these,  179  unique  data  were  measured  only  once  ,  829  twice  and 
31992  three  or  more  times.  We  note  that,  in  the  complete  set  of  unique  reflections,  25899 
reflections  were  found  to  be  greater  than  3ci(I).  The  internal  agreement  factors  based  on  the 
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averaging  of  symmetry  equivalent  reflections  are  R1  =  0.042,  R2  =  0.031,  vv^R  =  0.114  and 
S=1.4.  The  unit  cell  parameters,  data  collection,  averaging  and  other  crystallographic  parameters 
are  summarized  in  table  1 . 


3.  Conventional  Reflnement 

The  room  temperature  crystal  structure  of  16a,17P-estriol  has  been  previously  reported. 
We  solved  the  1 00  K  low  temperature  structure  by  direct  methods  using  the  program  SHELXS97 
and  the  spherical  atom  least  squares  refinement  was  carried  out  using  the  program  SHELXTL. 
All  the  hydrogen  atoms  were  located  from  a  difference  fourier  synthesis.  In  the  least  squares 
refinement,  all  non  hydrogen  atoms  were  treated  anisotropically  and  the  hydrogen  atoms  were 
isotropic.  The  refinement  converged  with  all  data  to  R=0.058,  wR=0.151  GOF=1.01. 
Comparison  of  Fq  and  Fc  values  of  the  largest  structure  factors  gave  no  indication  of  the  effect  of 
extinction  in  the  data.  A  thermal  ellipsoid  diagram  (figure  1 .)  shows  the  two  molecules  in  the 
asymmetric  unit  and  the  standard  atom  numbering  scheme. 

4.  Multipole  refinement 

In  order  to  model  the  deformation  electron  density  and  explore  the  bond  topological 
characteristics  like  p(r),  V^p(r),  bond  ellipticity  (s),  molecular  electrostatic  potential  etc.,  a 
multipole  aspherical  atom  refinement  was  carried  out  using  the  Hansen-Coppens  multipole 
formalism  implemented  in  the  XD  program.  According  to  this  formalism  the  total  electron 
density  of  an  aspherical  atom  is  parameterized  using  an  atom-centered  multipole  expansion 

patom(r)  =  Pc  pc  +  Pv  Pv(Kr)  +  E  k'/^  R/  (k'/F)!  ?im±  dim±  (0,(p) 


The  first  two  terms  pc  and  Pv(Kr)  are  the  spherically  averaged  Hartree-Fock  atomic  core  and 
valence  densities  normalized  for  one  electron.  The  third  term  describes  the  aspherical  part  of  the 
pseudo  atom.  R/(r)  =  [a/ "'^^/(ni+2)!]r'’'  e'®'"^  are  the  normalized  radial  distribution  functions,  a/ 
are  the  Hartree  Fock  optimized  single  zeta  values,  k  and  k’  are  radial  expansion  /contraction 
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parameters  whieh  are  optimized  in  the  least-square  refinement  along  with  Py  and  P/m-  P/m  are  the 
multipole  population  coefficients.  d/,„+  are  density  normalized  real  spherical  harmonic  angular 
functions  and  r,  0  and  cp  are  user  defined  local  atom  centered  coordinates  expressed  in  a  polar 
coordinate  system.  In  the  multipole  refinement  all  carbon  and  oxygen  atoms  were  allowed  to 
expand  up  to  the  hexadecapole  and  hydrogen  atoms  to  the  dipole  level.  The  k  and  k'  parameters 
of  all  non-H  atoms  were  refined,  but  hydrogen  atoms  k’  was  fixed  at  1.2. 

In  order  to  obtain  accurate  positional  and  thermal  parameters  for  the  non-hydrogen 
atoms,  a  high  angle  refinement  was  carried  out  using  the  data  with  sin0/X  >  0.7  A'*.  All  sp^,  sp^ 
and  hydroxyl  hydrogen  atoms  were  then  reset  to  neutron  bond  lengths  [Csp2-H=  1.083, 
Csp3-H=1.099(CH),  1.092(CH2),  1.059(CH3)  &  0-H  =  0.967A]  as  given  in  the  international 
table.  This  refinement  converged  without  any  significant  change  of  scale  factor,  thus  indicating 
the  quality  of  the  high  angle  data.  The  multipole  refinement  was  carried  out  on  25520  reflections 
[I>3a(I)]  with  the  resolution  (sin0/X,)max=l-lA’'.  The  refinement  strategy  was  as  follows:  In  the 
first  step  Py,  K  and  scale  factor  were  refined;  second  step  P/m,  k';  third  step  xyz,  Uy,  scale  factor 
were.  Finally  xyz,  Uy  and  all  electronic  parameters  were  refined  together.  Since  there  are  two 
molecules  in  the  asymmetric  unit,  their  multipole  parameters  were  constrained  to  be  equal  at  the 
beginning  of  the  refinement  except  the  0-H  group  atoms  as  these  are  involved  in  strong 
hydrogen  bonding.  In  addition,  to  avoid  the  phase  problem,  all  chemically  equivalent  atoms  were 
initially  constrained  to  be  equal  but  no  site  symmetry  was  imposed  on  any  atoms  in  the 
molecules.  At  the  later  stages  of  the  refinement,  the  chemieal  atom  constraints  within  the 
molecules  were  released  but  not  between  the  molecules.  The  residual  index  from  the  final 
refinement  is  given  in  table  1 .  An  attempt  was  made  to  refine  the  two  molecules  independently 
after  the  final  refinement.  This  refinement  did  not  give  a  meaningful  result  so  we  did  not  proceed 
further.  The  net  atomic  charges  were  calculated  from  a  K-refinement  (  Py,  k,  and  scale).  During 
the  refinements  the  unit  cell  was  constrained  to  be  neutral  and  no  charge  transfer  was  allowed 
between  the  molecules.  The  featureless  residual  density  maps  of  the  two  molecules  (figure  2) 
from  the  Fourier  summation  (Fo-Fmui)  indicates  the  adequacy  of  the  model  refinement. 

It  is  essential  to  test  the  physical  significance  and  the  correctness  of  the  thermal 
parameters  obtained  from  the  model  refinement.  To  accomplish  this,  Hirshfeld's  "rigid  bond" 
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test  was  performed  for  both  the  molecules  leading  to  a  satisfactory  result.  The  maximiun 
difference  Aa,b  between  the  atomic  mean  square  displacement  amplitudes  along  the  interatomic 
vectors  was  0.0012  A.  The  fractional  atomic  coordinates  and  thermal  parameters  are  given  in 
table  2  and  3. 


Result  and  Discussion 


1.  Structural  aspects 

The  report  of  the  room  temperature  crystal  structure  indicated  significant  conformational 
differences  between  the  two  molecules  in  the  asymmetric  unit.  We  observe  some  notable 
differences  between  the  reported  room  temperature  and  current  low  temperature  structure,  so  it  is 
worth  revisiting  the  structural  analysis.  Figure  1 .  shows  the  thermal  ellipsoid  representation  of 
estriol  with  two  independent  molecules  in  the  asymmetric  unit  obtained  from  the  low 
temperature  100  K  measurement.  The  geometrical  comparison  of  the  two  molecules  explicitly 
indicates  how  the  two  molecules  are  conformationally  different.  On  comparing  the  bond 
distances  between  the  two  molecules  no  significant  difference  were  observed  in  contrast  to  the 
previous  report.  The  maximum  deviation  in  the  bond  lengths  of  the  ring  (A-D)  systems  between 
the  two  molecules  is  0.01  A  (table  4a).  There  is  no  evidence  for  C(2)-C(3)  bond  shortening  in 
the  A-ring  systems  as  the  room  temperature  structure  reports.  The  average  Car-Car  bond 
distances  of  A-ring  atoms  of  molecule  1  and  2  are  1.3995  and  1.3990  A  respectively.  The  CH3 
group  of  the  second  molecule  lies  above  the  aromatic  ring  of  the  first  producing  distortion  in  the 
aromatic  ring.  This  effect  is  reproduced  here  producing  a  folding  in  the  A-ring  of  molecule  1 
with  a  dihedral  angle  of  2.99°.  There  is  no  such  effect  in  the  seeond  moleeule  where  the  A-ring 
system  is  planar  (0.52°).  There  are  differenees  in  bond  angles  between  the  molecules  1  and  2 
around  the  atoms  C(9),  C(13)  and  C(14)  [table  4b]  .  However,  we  observe  a  maximum  deviation 
of  2.5°  compared  to  3.9°  as  previously  reported.  The  torsion  angles  in  the  table  4  indicate  the 
distortion  of  the  B  and  C  ring  systems.  The  steric  hindrance  between  the  C-ring  equatorial 
hydrogen  atom  at  C(1 1)  and  the  hydrogen  atom  at  C(l)  produce  different  conformations  of  the  B 
and  C  ring  systems  in  the  two  molecules.  Comparison  of  torsion  angles  between  the  molecules 
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shows  the  major  distortion  appears  at  the  head  and  tail  part  of  the  moleeules  (table  5).  This  is 
attributed  to  the  different  environment  of  intermolecular  interactions  of  the  molecules  (1  and  2) 
with  their  neighbors.  The  relative  strengths  of  the  hydrogen  bonds  to  0(3)  in  the  two  molecules 
was  previously  discussed  on  the  basis  of  0.04  A  difference  in  the  C(16)-0(3)  bond  length.  We 
observe  less  than  0.01  A  difference  in  this  distance  between  molecules  1  and  2.  Hydrogen 
bonding  does,  however,  have  an  effect  on  the  geometry  of  the  hydroxyl  groups.  In  the  second 
molecule  the  16a  and  IVp-hydroxyl  C-O-H  bond  angles  are  systematically  less  compared  with 
the  first,  the  maximum  difference  being  8.8°.  There  is  less  of  a  differentiation  in  the  3-hydroxyl 
bonds,  the  deviation  being  2.1°.  The  structure  is  characterized  by  strong  hydrogen  bonding. 
Figure  2.  depicts  the  hydrogen  bonding  environment  of  the  two  molecules  in  the  asymmetric 
unit.  Each  hydroxyl  group  is  involved  in  at  least  two  hydrogen  bonds.  0(1’)  has  one  additional 
hydrogen  bond  resulting  in  a  tetrahedral  coordination  of  the  0-atom.  The  head  to  tail  hydrogen 
bonds  link  the  molecules  into  chains.  The  hydrogen  bond  geometrical  parameters  are  given  in 
table  6. 


2.  Charge  density 

The  model  static  deformation  density  map  in  figure  3  clearly  shows  the  bonding  features 
of  the  different  atoms  in  the  estriol  molecule.  All  chemical  bonds  show  an  increase  in  electron 
density  with  respect  to  neutral  spherical  atoms  as  well  as  density  associated  with  the  lone  pair 
regions  of  the  oxygen  atoms.  The  featureless  residual  electron  density  map  was  computed  using 
the  XDFOUR  program  after  the  multipole  refinement.  The  featureless  residual  maps  (figure  2a 
and  b)  of  the  two  molecules  in  the  asymmetric  unit,  confirms  the  validity  and  the  adequacy  of 
the  model  refinement  and  the  quality  of  the  data  set.  The  minimum  and  maximum  residual 
density  calculated  for  the  molecules  are  0.21,  -0.29  (moleculel)  and  0.22,  -0.27eA'^  (molecule 
2)  respectively.  This  values  are  very  close  to  the  requirement  of  the  average  error  (±0.05)  eA'^  in 
the  experimental  difference  density  maps.  As  specified  earlier,  in  the  multipole  refinement,  the 
respective  chemically  equivalent  ring  atoms  of  the  two  molecules  were  constrained  pair  wise. 
No  constraint  was  imposed  on  the  hydroxyl  groups,  because  they  involved  strong  hydrogen 
bonding.  All  further  discussions  are  based  on  this  constraint  model.  In  order  to  further 
characterize  the  interaction  of  atoms  in  the  molecule  a  bond  topological  analysis  was  performed 
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based  on  Bader’s  theory  of  “atoms  in  molecules”  implemented  in  XDPROP.  This  analysis 
includes  the  determination  of  electron  density  pc(r)  at  the  bond  critical  point  (CP)  where  Vp(r)  = 
0  between  the  bonded  atoms,  its  Laplacian  V^pc(r)  and  bond  ellipticity  [s  =  (X,i/X2)-1].  Also,  the 
eigenvalues  of  the  second  derivatives  of  p(r)  was  calculated  to  characterize  the  nature  of  the 
bond  critical  point.  The  bond  path  analysis  gives  information  about  the  amount  of  deviation  of 
the  bond  critical  point  from  the  intemuclear  axis,  thus  obtaining  an  additional  measure  of  strain. 
All  these  results  including  the  length  of  bond  paths  and  intemuclear  distances  are  listed  in  table 
7. 

A-ring  system:  The  topological  analysis  results  reported  in  table  7,  such  as  A,i_  A,2,  ^,3,  pc(r)  , 
V^pc(r),  8  reveal  electronic  insight  into  the  estriol  molecules.  The  magnitude  of  pc(r)  at  the  bond 
critical  point  is  a  clear  indicator  of  chemical  bond  strength.  The  electron  density  pc(r)  at  the 
critical  points  of  the  aromatic  ring  (Csp2  -Csp2 )  bonds  range  from  2.04  to  2.16  eA’^.  The  average 
value  (2.1 1  eA'^)  is  a  typical  value  for  aromatic  bonds.  The  Laplacian  V^pc(r)  is  a  very  sensitive 
parameter  which  identifies  the  areas  where  charge  is  concentrated  (negative  value)  or  depleted 
(positive  value)  in  the  bonding  region.  In  this  respect  we  observe  the  maximum  negative 
Laplacian  in  the  A-ring  bonds  is  -22.2(1)  eA'^  and  this  value  appears  less  in  the  bonds  linked  to 
B  ring  atoms.  However  the  average  value  is  -18.72  eA'^  which  is  very  close  to  an  earlier 
observation  for  the  benzene  molecule.  The  dispersion  of  bond  ellipticity  is  attributed  to  the 
preferential  charge  accumulation  in  the  bond.  Bader’s  theory  on  ellipticity  emphasis  the  concept 
of  the  a-rr  character  of  the  bond.  Explicitly,  it  is  the  measure  of  deviation  of  charge  density  from 
spherical  symmetry  (at  the  critical  point).  In  this  aspect  the  average  value  (0.17)  of  ellipticity  of 
Csp2  -Csp2  bonds  indicates  significant  n  (aromatic)  character.  The  Pc(r)’s  in  the  Csp2  -H  bonds 
ranges  from  1.78(4)  to  1.89(1)  eA'^  and  this  density  is  consistent  with  one  would  expect  for 
single  C-H  bonds.  The  Laplacian  the  of  C(l)-H(l)  bond  is  small  compared  with  the  other  two  C- 
H  bonds,  perhaps  indicating  weakening  of  the  bond  due  to  steric  interactions  with  hydrogen 
atoms  at  C(1 1).  The  small  ellipticity  of  Csp2-H  bonds  is  consistent  with  only  a  ct  contribution  in 
the  single  bonds. 

B  and  C-  ring  systems:  The  density  of  chemically  equivalent  Csp2  -Csps  bonds  in  the  B-ring 
are  equal  and  their  Laplacian’ s  are  very  close  (see  the  table  7).  In  both  ring  systems,  the  average 
values  (  1.70  and  1.68  eA'^ )  of  Pc(r)’s  in  the  Csp3  -Csp3  bonds  are  not  very  different  and  a 
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similar  trend  appears  for  the  Laplacian.  The  pc(r)’s  and  Laplacian  for  these  bonds  are 
significantly  less  than  the  aromatic  ring  system  in  agreement  with  their  single  bond  character.  In 
the  C-ring  the  ellipticity  of  the  bonds  C(ll)-C(12)  and  C(12)-C(13)  are  unexpectedly  larger 
compared  to  the  similar  bonds  having  small  values  (see  the  table  7)  indeed  the  values  approach 
those  of  the  aromatic  ring.  In  all  Csp3  -H  (CH2)  bonds  the  density  ranges  from  1.64  to  1.87 
eA'^  and  Csp3  -H  (  CH)  bonds  1.75  to  1.77  eA'^  The  average  effect  is  same  in  both  types  of 
bonds.  The  Laplacian  of  these  bonds  is  consistently  smaller  than  for  the  Csp2  -H  bonds.  The 
maximum  and  minimum  values  are  -8.8(2),  -1 8.3(1)  eA’^ . 

D-ring  system:  Due  to  this  being  a  5  membered  ring  and  to  the  different  hydrogen  bonding 
environments  this  ring  is  highly  strained  in  the  geometry  of  both  molecules  (1  and  2).  Although 
its  geometry  is  strained,  no  significant  difference  was  observed  in  the  Pc(r)’s  at  the  bond  critical 
points  except  for  C(14)-C(15)  ==  1.57  eA’^  implying  a  weakening  of  this  bond.  The  average  value 
of  the  Laplacian  is  -12.78  eA'^  We  observe  the  ellipticity  of  this  five  membered  system  is  close 
to  B-ring  bonds  and  they  are  small  consistent  with  single  bonds.  The  electron  density  (1.95  eA'^) 
and  Laplacian  of  C(16)-H(16)  bond  is  surprisingly  higher  than  the  similar  C(17)-H(17)  bond  in 
this  ring  (see  the  table  7). 

C-O  and  O-H  groups:  These  groups  play  a  significant  role  in  the  whole  geometry  of  both  the 
molecules.  The  Pc(r)’s  of  Csp2  -  O  [A-rings]  bonds  of  both  the  molecules(  1  and  2)  are  larger 
than  Csp3  -O  bonds  of  D-ring  in  agreement  with  conjugation  of  0(1)  to  the  aromatic  ring.  The 
average  values  are  2.1 1  and  1.91  eA  ^ .  The  Laplacians  of  the  Csp2  ~  O  bonds  are  not  equal  to 
Csp3  -  O  bonds  and  also  indicate  a  strengthening  of  the  bond  to  the  aromatic  ring.  Whereas  the 
ellipticity  of  the  0(1)-C(3)  bond  in  molecule  1  is  a  strong  indicator  of  a  7t  component  to  the 
bond,  that  in  molecule  2  is  closer  to  the  value  observed  for  the  D  ring.  We  attribute  this  to  the 
difference  in  the  hydrogen  bonding  in  the  two  molecule. 

The  small  values  at  the  ellipticity  (0.06-0.1)  in  the  Csp3  -  O  bonds  for  both  molecules  indicates 
an  unperturbed  pure  ct  bonds. 

Hydrogen  bonding:  Table  8  reports  the  hydrogen  bond  topology.  There  are  three  intermolecular 
hydrogen  bonds  in  each  molecule  that  meet  Bader’s  criteria  for  polar  non-bonded  interaction  of 
the  type  H-  •  •  A  by  observing  a  CP  along  the  H-  •  -A  path.  Of  the  three  hydrogen  bonds  from 
molecule  1,  the  two  bonds  [0(1)-H(1A)-  •  •  0(3)  and  0(3)-H(3A)-  •  •  0(1’)]  display 
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approximately  linear  hydrogen  bonding,  171.7,  170.5°  respectively.  The  corresponding  electron 
densities  of  these  bonds  at  the  CP  are  similar  with  values  of  0.23  and  0.20  eA’^  respectively. 
Both  of  these  bonds  have  a  positive  Laplacian  V^p  at  the  CP.  This  is  characteristic  of  closed 
shell  interactions  which  are  governed  by  the  contraction  of  the  electro  density  towards  each  of 
the  interacting  atoms.  In  molecule  2  we  observe  similar  electron  densities  and  Laplacians  for  the 
bonds  0(r)-H(rA)-  •  •  0(3’)  and  0(2’)-H(2’A)-  •  •  0(2’).  The  hydrogen  bond  0(2’)-H(2’A)-  • 

•  0(2’)  in  molecule  1  shows  moderate  electron  density  at  the  CP  (0.11  eA‘^)  with  a  positive 
Laplacian.  Although  the  hydrogen  bond  0(2’)-H(2’A)-  •  •  0(1’)  is  approximately  linear  (169.2°) 
the  electron  density  is  very  small  at  the  CP  (0.05  eA'^).  This  may  be  attributed  to  the  longer 
H(2’A)-  •  -0(1  ’A)  distance  and  the  fact  that  O(rA)  accepts  one  additional  hydrogen  bond. 

3.  Electrostatic  Potential 

An  additional  important  application  of  the  multipole  model  is  that  it  provides  an 
opportunity  to  understand  the  chemical  reactivity  of  molecules.  That  is,  to  identify  the  places 
where  nucleophilic  or  electrophilic  attack  may  take  place.  The  regions  of  positive  potential  will 
attract  nucleophilic  reagents  and  the  electronegative  regions  of  the  molecule  will  determine  the 
approach  of  electrophilic  reagents.  The  electrostatic  potential  (ESP)  not  only  plays  an  important 
role  in  directing  chemical  attack,  but  also  for  molecular  recognition  ( drug  receptor  interaction). 
This  is  because  electrostatic  forces  are  long  range  interactions  which  determine  the  most  stable 
geometrical  orientation  of  two  molecules.  To  investigate  this  aspect  we  have  calculated  the 
electrostatic  potential  around  the  molecules  by  the  method  of  Su  and  Coppens  based  on  the 
populations  obtained  from  the  multipole  refinement.  The  figures  4a  and  4b  show  a  contour  plot 
of  the  eleetrostatic  potential  of  the  molecule  1  and  2  in  the  molecular  plane.  The  positive 
contours  around  the  aromatic  rings  are  more  diffuse  than  those  of  the  saturated  rings.  The 
negative  contours  around  the  16a,  17p  oxygen  atoms  are  much  more  pronounced  than  for  the  A 
ring.  The  potential  around  molecule  1  is  more  symmetric  than  for  molecule  2.  Figure  5  shows  the 
electrostatic  potential  in  the  hydrogen  bonding  environment  of  each  oxygen  atom  of  the 
molecule. 

Figure  6  displays  the  3D  isosurface  representation  of  the  electrostatic  potential  computed 
from  the  molecule  extracted  from  the  crystal.  The  strong  extended  electro  negative  regions  are 
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evident  at  the  head  and  tail  portion  of  the  molecule.  The  shape  and  direction  of  polarization  is 
attributed  to  the  intermolecular  interaction. 

4.  Net  atomic  charges 

The  net  atomic  charges  (table  9)  were  calculated  from  from  a  K-refinement  (  k,  Pv,  scale). 
In  this  refinement  the  chemically  equivalent  atoms  were  constrained  as  described  in  the 
multipole  refinement.  All  oxygen  atoms  carry  significant  negative  charges  in  the  range  -0.63  to  - 
0.87e.  Similarly  all  hydroxyl  hydrogen  atoms  carry  high  positive  charges,  approximately  double 
the  amount  of  charge  of  the  other  hydrogen  atoms  in  the  molecules.  A  group  charge  calculation 
on  both  molecules  indicates  that  the  central  regions  are  close  to  neutral,  whereas  the  head  and  tail 
regions  form  charged  units. 


Conclusion 

The  present  study  on  estriol  gives  information  on  both  the  structural  and  electronic 
properties  of  the  molecule.  The  differences  in  hydrogen  bonding,  molecular  packing 
environment,  and  interamolecular  steric  effects  are  responsible  for  small  differences  in  the  two 
molecules  in  the  asymmetric  unit.  The  significant  differences  as  reported  for  the  room 
temperature  structure  were  not  confirmed  here.  The  comparison  of  torsion  angles  between  the 
two  molecules  shows  that  the  major  differences  at  the  head  and  tail  portions  of  the  molecules 
may  be  attributed  to  the  different  hydrogen  bonding  environment,  both  molecules  being  involved 
in  strong  hydrogen  bonding. 

The  electron  densities  at  the  CP  of  the  aromatic  C-C  bonds  is  larger  than  for  the  Csps 
-Csp3  bonds  of  the  B,C  and  D  ring  systems.  The  additional  electron  density  accumulation  in  the 
aromatic  bonds  along  with  the  ellipticity  indicates  a  n  contribution.  The  unconstrained 
refinement  of  the  hydroxyl  groups  between  the  two  molecules  shows  differences  in  both  C-0 
and  0-H  bond  electron  densities  at  the  CP.  The  average  Csp3  -O  bond  densities  significantly 
less  than  Csp2-0.  The  charge  accumulation  in  the  0-H  (16a,17p  position)  bond  is  slightly  larger 
than  in  the  3-hydroxy  0-H  bond.  The  bond  path  calculation  reveals  a  deviation  of  the  CP  from 
the  intemuclear  axis  for  several  C-0  bonds.  The  CP  of  Csp2-0  bond  in  the  first  molecule  deviates 
by  0.03  A  from  the  intemuclear  axis,  while  in  molecule  2  it  lies  very  close  (0.006  A).  In 
molecule  1  C(17)-0(2)  and  C(16)-  0(3)  deviate  0.03  and  O.OlA  respectively.  But  these 
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deviations  are  exactly  opposite  in  the  bonds  C(17’)-0(2’)  [  O.OlA]  and  C(16’)-  0(3’)  [0.03A]  of 
molecule  2.  Strong  electronegative  regions  were  found  in  the  electrostatic  potential  maps  at  the 
vicinity  of  the  oxygen  atoms  of  both  molecules.  The  extension  of  positive  and  negative  contours 
in  both  molecules  explicitly  indicates  the  electrostatic  interaction  with  the  neighboring 
molecules.  The  large  electron  density  at  the  hydrogen  bond  CP’s  implies  the  very  strong  nature 
of  hydrogen  bonding  in  the  estriol  crystal.  Net  atomic  charges  computed  from  the  K-refinement 
gave  a  reasonable  charge  distribution.  The  oxygen  atoms  0(2), 0(3)  (molecule  1)  and 
0(2’), 0(3’)  (molecule  2)  carry  a  higher  negative  charge  than  0(1)  and  0(1  ’). 
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Table  1.  Crystal  data  and  experimental  conditions 

Crystal  data 

Empirical  formula 

C18H24O3 

Formula  weight 

288.37 

Crystal  system 

Monoclinic 

Space  group 

P2, 

a  (A) 

7.5077(4) 

b(A) 

23.0809(12) 

c  (A) 

9.1632(5) 

PO 

110.957(2) 

V(A') 

1482.8 

z 

4 

F(OOO) 

624 

Dx(Mg/m^) 

1.292 

p  (mm'') 

0.055 

Data  collection 

Diffractometer 

Bruker’s  SMART  2K  CCD 

Radiation,  X  (A) 

Ag  Ka,  0.56087 

Scan  method 

o-scan 

No.  of  reflections  measured 

278184 

Sin(e/A.)^ax  (A-') 

1.29 

Rint 

0.042 

No.  of  unique  data 

43918 

Refinement 

No.  of  relections  (F>3ct,  Sin(0/X)max=l-1  A'') 

25520 

R  (F),  wR(F) 

0.031,0.035 

GOF 

0.98 

Nobs/Npar 

27.3 

35 


Table  2.  Fractional  coordinates,  equivalent  displacement  parameters  after 
multipole  refinement. 


Atom 

X 

y 

z 

Ue, 

0(1) 

0.51770(10) 

0.85356(11) 

0.52121(11) 

0.018 

0(2) 

0.79922(9) 

1.32165(11) 

0.40946(9) 

0.015 

0(3) 

1.19063(9) 

1.28093(11) 

0.40370(9) 

0.015 

0(1') 

0.98786(9) 

1.26564(11) 

1.07236(8) 

0.013 

0(2') 

0.74949(10) 

0.79559(1 1) 

1.08082(9) 

0.016 

0(3') 

0.31380(9) 

0.81946(11) 

0.90942(8) 

0.015 

C(l) 

0.45781(8) 

1.00734(11) 

0.42503(8) 

0.014 

C(2) 

0.41441(8) 

0.94973(11) 

0.44369(8) 

0.016 

C(3) 

0.56256(8) 

0.91063(11) 

0.51232(7) 

0.012 

C(4) 

0.75074(7) 

0.93008(11) 

0.56559(6) 

0.011 

C(5) 

0.79387(7) 

0.98793(11) 

0.54672(6) 

0.010 

C(6) 

1.00171(8) 

1.00583(11) 

0.60532(7) 

0.013 

C(7) 

1.03656(7) 

1.06196(11) 

0.53177(7) 

0.013 

C(8) 

0.89490(7) 

1.10848(11) 

0.53903(6) 

0.010 

C(9) 

0.69129(7) 

1.08934(11) 

0.43668(6) 

0.010 

C(10) 

0.64663(7) 

1.02776(11) 

0.47412(6) 

0.010 

C(ll) 

0.54264(8) 

1.13511(11) 

0.43758(7) 

0.013 

C(12) 

0.58776(7) 

1.19498(11) 

0.38487(7) 

0.012 

C(13) 

0.78968(7) 

1.21489(11) 

0.48383(6) 

0.010 

C(14) 

0.93243(7) 

1.16705(11) 

0.47887(6) 

0.010 

C(15) 

1.12818(8) 

1.19669(11) 

0.55148(8) 

0.015 

C(16) 

1.08802(7) 

1.26095(11) 

0.49991(7) 

0.012 

C(17) 

0.87053(7) 

1.26406(11) 

0.41297(6) 

0.010 

C(18) 

0.80161(10) 

1.23053(11) 

0.65010(7) 

0.017 

0(1') 

1.04605(7) 

1.11001(11) 

1.03189(7) 

0.013 

0(2') 

1.08893(7) 

1.16882(11) 

1.05868(7) 

0.013 

0(3') 

0.94251(7) 

1.20772(11) 

1.04701(6) 

0.010 

0(4') 

0.75671(7) 

1.18775(11) 

1.01090(6) 

0.010 

0(5') 

0.71459(7) 

1.12879(11) 

0.98527(6) 

0.009 

0(6') 

0.51022(7) 

1.11048(11) 

0.95314(8) 

0.013 

0(7') 

0.47065(7) 

1.04732(11) 

0.90308(7) 

0.012 

0(8') 

0.63561(7) 

1.00870(11) 

0.99990(6) 

0.009 

0(9') 

0.81329(7) 

1.02458(11) 

0.96111(6) 

0.009 

C(IO') 

0.86000(7) 

1.08864(11) 

0.99458(6) 

0.010 

0(1 1') 

0.98225(7) 

0.98344(11) 

1.03926(7) 

0.012 

C(12') 

0.92731(7) 

0.91897(11) 

1.00723(7) 

0.012 

C(13') 

0.75706(7) 

0.90430(11) 

1.05525(6) 

0.009 

C(14') 

0.59062(7) 

0.94491(11) 

0.96631(6) 

0.009 

0(15') 

0.41988(8) 

0.91953(11) 

0.99962(8) 

0.014 

0(16') 

0.46121(8) 

0.85359(11) 

1.01718(7) 

0.012 

C(17') 

0.65667(8) 

0.84561(11) 

0.99893(6) 

0.011 

C(18') 

0.81623(9) 

0.90818(11) 

1.23339(6) 

0.015 

H(1A) 

0.63209 

0.82990 

0.55330 

0.034 

H(2A) 

0.69428 

1.32373 

0.44686 

0.012 

H(3A) 

1.13395 

1.27295 

0.29273 

0.025 

H(l) 

0.34450 

1.03712 

0.36414 

0.029 

H(2) 

0.26604 

0.93672 

0.39812 

0.044 

H(4) 

0.86510 

0.89941 

0.61941 

0.021 

H(6A) 

1.09283 

0.97171 

0.58700 

0.037 
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H(6B) 

1.04093 

1.01387 

0.73024 

0.020 

H(7A) 

1.01476 

1.05396 

0.40896 

0.015 

H(7B) 

1.18207 

1.07730 

0.59177 

0.021 

H(8) 

0.89840 

1.11151 

0.65964 

0.017 

H(9) 

0.68554 

1.08694 

0.31532 

0.016 

H(llA) 

0.40304 

1.12061 

0.35796 

0.024 

H(llB) 

0.53292 

1.13757 

0.55336 

0.020 

H(12A) 

0.57626 

1.19302 

0.26278 

0.018 

H(12B) 

0.48167 

1.22649 

0.38902 

0.019 

H(14) 

0.90581 

1.15747 

0.35520 

0.014 

H(15A) 

1.23477 

1.17769 

0.51002 

0.022 

H(I5B) 

1.18913 

1.19417 

0.67908 

0.024 

H(16) 

1.12808 

1.29047 

0.60016 

0.025 

H{17) 

0.84636 

1.25096 

0.29187 

0.011 

H(18A) 

0.76328 

1.19583 

0.70939 

0.033 

H(18B) 

0.70167 

1.26365 

0.64136 

0.026 

H(18C) 

0.94099 

1.24472 

0.71651 

0.025 

H(rA) 

0.87426 

1 .28653 

1.06698 

0.027 

H(2'A) 

0.84956 

0.78829 

1.03968 

0.022 

H(3'A) 

0.28464 

0.83420 

0.80471 

0.020 

H(l') 

1.15901 

1.07987 

1.03660 

0.026 

H(2') 

1.23032 

1.18590 

1.08056 

0.026 

H(4’) 

0.64298 

1.21784 

1.00111 

0.021 

H(6'C) 

0.41147 

1.14015 

0.87102 

0.027 

H(6’D) 

0.48205 

1.11487 

1.06167 

0.014 

H(7'C) 

0.44868 

1.04194 

0.77938 

0.024 

H(7'D) 

0.33470 

1.03427 

0.91168 

0.021 

H(8') 

0.66573 

1.01546 

1.12513 

0.012 

H(9’) 

0.77340 

1.01922 

0.83423 

0.011 

H(llC) 

1.09205 

0.99035 

0.98873 

0.020 

H(nD) 

1.04447 

0.99134 

1.16493 

0.017 

H(12C) 

0.89389 

0.91066 

0.88280 

0.018 

H(12D) 

1.05438 

0.89402 

1.07392 

0.029 

H(14') 

0.56308 

0.93793 

0.84142 

0.012 

H(15C) 

0.28167 

0.92827 

0.90891 

0.023 

H(15D) 

0.40421 

0.93935 

1.10220 

0.016 

H(16') 

0.48007 

0.83935 

1.13620 

0.015 

H(17') 

0.62869 

0.84147 

0.87314 

0.014 

H(18D) 

0.86807 

0.94995 

1.27305 

0.017 

H(18E) 

0.92207 

0.87636 

1.28493 

0.022 

H(18F) 

0.70244 

0.89942 

1.27329 

0.032 

Table  3.  Anisotropic  thermal  parameters  of  non-Hydrogen  atoms  (A^)  and  estimated 
standard  deviations. 


Atom  U||  U22  tJ33  U12  U|3  U23 


0(1)  0.0125(2)  0.0088(2) 

0(2)  0.0146(2)  0.0087(2) 

0(3)  0.0113(2)  0.0137(2) 

0(1')  0.0118(2)  0.0071(1) 

0(2')  0.0171(2)  0.0093(2) 

0(3')  0.0144(2)  0.0111(2) 

C(l)  0.0083(2)  0.0112(2) 
C(2)  0.0094(2)  0.0109(2) 

C(3)  0.0108(2)  0.0089(1) 

C(4)  0.0103(2)  0.0096(1) 

C(5)  0.0087(1)  0.0086(1) 

C(6)  0.0085(2)  0.0104(2) 

C(7)  0.0083(2)  0.0104(2) 

C(8)  0.0080(1)  0.0089(1) 

C(9)  0.0077(1)  0.0086(1) 

C(10)  0.0080(1)  0.0090(1) 
C(ll)  0.0085(2)  0.0099(2) 
C(12)  0.0086(1)  0.0094(2) 
C(13)  0.0091(1)  0.0085(1) 
C(14)  0.0078(1)  0.0087(1) 
C(15)  0.0082(2)  0.0115(2) 
C(16)  0.0097(2)  0.0096(2) 
C(17)  0.0097(2)  0.0080(1) 
C(18)  0.0216(2)  0.0145(2) 
C(r)  0.0082(1)  0.0087(1) 
C(2')  0.0092(2)  0.0083(2) 

C(3')  0.0097(1)  0.0073(1) 

C(4')  0.0089(1)  0.0069(1) 

C(5')  0.0074(1)  0.0070(1) 

C(6')  0.0071(1)  0.0089(1) 

C(7')  0.0075(1)  0.0087(2) 

C(8')  0.0072(1)  0.0073(1) 

C(9')  0.0085(1)  0.0068(1) 

C(IO')  0.0078(1)  0.0076(1) 
C(ir)  0.0080(1)  0.0085(1) 
C(12')  0.0103(2)  0.0083(1) 
C(13')  0.0090(1)  0.0083(1) 
C(14')  0.0080(1)  0.0073(1) 
C(15')  0.0099(2)  0.0091(2) 
C(16')  0.0119(2)  0.0080(1) 
C(17')  0.0117(2)  0.0076(1) 
C(18')  0.0175(2)  0.0149(2) 


0.0312(3) 

0.0217(3) 

0.0187(2) 

0.0182(2) 

0.0200(3) 

0.0168(2) 

0.0211(2) 

0.0250(3) 

0.0164(2) 

0.0120(2) 

0.0103(2) 

0.0167(2) 

0.0170(2) 

0.0104(2) 

0.0108(2) 

0.0112(2) 

0.0193(2) 

0.0162(2) 

0.0106(2) 

0.0110(2) 

0.0214(2) 

0.0147(2) 

0.0123(2) 

0.0123(2) 

0.0217(2) 

0.0211(2) 

0.0128(2) 

0.0134(2) 

0.0125(2) 

0.0212(2) 

0.0175(2) 

0.0115(2) 

0.0112(2) 

0.0123(2) 

0.0174(2) 

0.0156(2) 

0.0089(1) 

0.0113(2) 

0.0220(2) 

0.0144(2) 

0.0121(2) 

0.0095(2) 


-0.0001(2) 

0.0028(1) 

-0.0035(2) 

-0.0011(1) 

0.0028(2) 

-0.0051(2) 

0.0005(1) 

-0.0003(1) 

0.0002(1) 

0.0005(1) 

0.0005(1) 

0.0010(1) 

0.0018(1) 

0.0008(1) 

0.0011(1) 

0.0005(1) 

0.0014(1) 

0.0014(1) 

0.0005(1) 

0.0011(1) 

0.0006(1) 

-0.0005(1) 

0.0010(1) 

-0.0001(2) 

0.0000(1) 

-0.0007(1) 

-0.0004(1) 

0.0002(1) 

0.0004(1) 

0.0007(1) 

0.0003(1) 

-0.0001(1) 

0.0003(1) 

0.0004(1) 

0.0008(1) 

0.0014(1) 

0.0000(1) 

-0.0003(1) 

-0.0015(1) 

-0.0017(1) 

-0.0003(1) 

-0.0022(2) 


0.0022(2) 

0.0011(2) 

-0.0017(2) 

-0.0007(1) 

0.0024(2) 

-0.0028(2) 

0.0020(2) 

0.0016(2) 

0.0008(1) 

0.0009(1) 

0.0007(1) 

0.0032(1) 

0.0024(1) 

0.0010(1) 

0.0008(1) 

0.0007(1) 

0.0017(1) 

0.0016(1) 

-0.0005(1) 

0.0008(1) 

0.0023(2) 

-0.0009(1) 

0.0001(1) 

-0.0018(2) 

-0.0009(2) 

-0.0005(1) 

-0.0003(1) 

0.0002(1) 

0.0001(1) 

-0.0006(1) 

-0.0019(1) 

-0.0016(1) 

-0.0006(1) 

0.0000(1) 

0.0005(1) 

0.0002(1) 

-0.0007(1) 

-0.0021(1) 

-0.0035(2) 

-0.0019(1) 

-0.0013(1) 

-0.0009(1) 


0.0108(2) 

0.0102(2) 

0.0066(2) 

0.0073(2) 

0.0096(2) 

0.0065(2) 

0.0039(2) 

0.0062(2) 

0.0071(1) 

0.0048(1) 

0.0029(1) 

0.0015(1) 

0.0041(1) 

0.0024(1) 

0.0023(1) 

0.0026(1) 

0.0054(1) 

0.0034(1) 

0.0042(1) 

0.0028(1) 

0.0020(2) 

0.0041(1) 

0.0048(1) 

0.0089(2) 

0.0066(1) 

0.0067(2) 

0.0053(1) 

0.0041(1) 

0.0032(1) 

0.0034(1) 

0.0008(1) 

0.0034(1) 

0.0041(1) 

0.0040(1) 

0.0047(1) 

0.0069(1) 

0.0040(1) 

0.0032(1) 

0.0077(2) 

0.0066(1) 

0.0060(1) 

0.0044(2) 
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Table  4a.  Selected  bond  lengths  (A) 


Atoms 

C(l)-C(2) 

C(2)-C(3) 

C(3)-C(4) 

C(4)-C(5) 

C(5)-C(10) 

C(l)-C(10) 

0(1)-C(3) 

C(5)-C(6) 

C(6)-C(7) 

C(7)-C(8) 

C(8)-C(9) 

C(9)-C(10) 

C(9)-C(ll) 

C(ll)-C(12) 

C(12)-C(13) 

C(13)-C(14) 

C(13)-C(17) 

C(13)-C(18) 

C(14)-C(15) 

C(14)-C(8) 

C(15)-C(16) 

C(16H17) 

C(16)-0(3) 

C(17)-0(2) 


Molecule(l) 

1.3939(8) 

1.3953(8) 

1.3938(7) 

1.3991(7) 

1.4084(7) 

1.4063(7) 

1.3692(8) 

1.5147(7) 

1.5256(7) 

1.5293(7) 

1.5446(7) 

1.5274(7) 

1.5389(7) 

1.5408(7) 

1.5324(7) 

1.5508(7) 

1.5360(7) 

1.5368(8) 

1.5399(7) 

1.5237(7) 

1.5529(7) 

1.5412(7) 

1.4386(9) 

1.4291(7) 


Molecule(2) 

1.3962(7) 

1.3936(7) 

1.3920(7) 

1.3978(6) 

1.4110(7) 

1.4035(7) 

1.3784(7) 

1.5154(7) 

1.5252(7) 

1.5263(7) 

1.5425(6) 

1.5252(6) 

1.5407(7) 

1.5442(7) 

1.5311(7) 

1.5413(7) 

1.5464(7) 

1.5333(7) 

1.5360(7) 

1.5172(7) 

1.5501(7) 

1.5462(7) 

1.4278(8) 

1.4167(8) 


Table  4b.  Selected  bond  angles  (°) 


Bonds 

Molecule  (1)  Molecule(2) 

C(2)- 

0(1)- 

0(10) 

122.2(1) 

122.1(1) 

C(l)- 

0(2)- 

0(3) 

119.2(1) 

199.1(1) 

0(1)- 

0(3)- 

0(2) 

118.4(1) 

118.1(1) 

0(1)- 

0(3)- 

0(4) 

122.0(1) 

121.8(1) 

C(2)- 

0(3)- 

0(4) 

119.6(1) 

120.0(1) 

C(3)- 

0(4)- 

0(5) 

121.0(1) 

120.6(1) 

C(4)- 

0(5)- 

0(6) 

118.1(1) 

117.2(1) 

C(4)- 

0(5)- 

0(10) 

120.3(1) 

120.4(1) 

C(6)- 

0(5)- 

0(10) 

121.6(1) 

122.4(1) 

C(5)- 

0(6)- 

0(7) 

113.1(1) 

113.6(1) 

C(6)- 

0(7)- 

0(8) 

110.2(1) 

110.7(1) 

C(7)- 

0(8)- 

0(9) 

109.0(1) 

108.4(1) 

C(7)- 

0(8)- 

0(14) 

112.7(1) 

111.9(1) 

C(9)- 

0(8)- 

0(14) 

107.8(1) 

109.4(1) 

C(8)- 

0(9)- 

0(10) 

112.3(1) 

109.9(1) 

C(8)- 

0(9)- 

0(11) 

111.3(1) 

112.5(1) 

C(10) 

-  C(9)  - 

0(11) 

114.8(1) 

114.3(1) 

C(l)- 

0(10)- 

-0(5) 

117.6(1) 

117.7(1) 

C(l)- 

0(10). 

-0(9) 

121.1(1) 

121.8(1) 

C(5)- 

0(10). 

-0(9) 

121.0(1) 

120.5(1) 

C(9)- 

0(11)- 

-0(12) 

111.9(1) 

112.7(1) 

C(ll) 

-0(12) 

-0(13) 

111.4(1) 

110.6(1) 

C(12) 

-0(13) 

-0(14) 

108.5(1) 

108.4(1) 

C(12) 

-0(13) 

-0(17) 

115.4(1) 

117.1(1) 

C(12) 

-0(13) 

-0(18) 

110.6(1) 

109.8(1) 

0(14) 

-0(13) 

-0(17) 

97.6(1) 

98.6(1) 

0(14) 

-0(13) 

-0(18) 

113.6(1) 

113.5(1) 

0(17) 

-0(13) 

-0(18) 

110.7(1) 

109.1(1) 

0(8)- 

0(14). 

-0(13) 

113.7(1) 

113.8(1) 

0(8)- 

0(14) . 

-0(15) 

120.6(1) 

118.1(1) 

0(13) 

-  C(14) 

-0(15) 

103.5(1) 

103.9(1) 

0(14) 

-0(15) 

-0(16) 

104.0(1) 

104.2(1) 

0(3)- 

0(16) 

-0(15) 

113.9(1) 

113.1(1) 

0(3)- 

0(16) 

-0(17) 

111.6(1) 

113.0(1) 

0(15) 

-0(16) 

-0(17) 

104.8(1) 

105.6(1) 

0(2)- 

0(17) 

-0(13) 

119.7(1) 

116.6(1) 

0(2)- 

0(17) 

-0(16) 

112.0(1) 

110.2(1) 

0(13) 

-  C(16) 

-0(17) 

104.6(1) 

103.2(1) 

0(3)- 

0(1)- 

H(1A) 

110.5(2) 

108.4(2) 

0(17) 

-  0(2)  - 

H(2A) 

112.6(2) 

103.8(2) 

0(16) 

-  0(3)  ■ 

•  H(3A) 

116.8(1) 

109.7(2) 

40 


Table  4b  (continued).  H-atom  bond  angles  (°) 


Bonds  Molecule  (1)  Molecule  (2) 


C(3)-  0(1)-  H(1A) 
C(17)  -  0(2)  -  H(2A) 
C(16)  -  0(3)  -  H(3A) 
C(2)-  C(l)-  H(l) 
C(10)-C(l)-  H(l) 
C(l)-  C(2)-  H(2) 

C(3)  -  C(2)  -  H(2) 

C(3)-  C(4)-  H(4) 
C(5)-  C(4)-  H(4) 

C(5)  -  C(6)  -  H(6A) 
C(5)  -  C(6)  -  H(6B) 
C(7)  -  C(6)  -  H(6A) 
C(7)  -  C(6)  -  H(6B) 
H(6A)  -  C(6)  -  H(6B) 
C(6)  -  C(7)  -  H(7A) 
C(6)  -  C(7)  -  H(7B) 
C(8)  -  C(7)  -  H(7A) 
C(8)  -  C(7)  -  H(7B) 
H(7A)  -  C(7)  -  H(7B) 
C(7)  -  C(8)  -  H(8) 

C(9)  -  C(8)  -  H(8) 
C(14)  -  C(8)  -  H(8) 

C(8)  -  C(9)  -  H(9) 

C(10)-C(9)-  H(9) 

C(ll)-C(9)-  H(9) 

C(9)-  C(11)-H(11A) 

C(9)-  C(11)-H(11B) 

C(12)-C(11)-H(11A) 

C(12)-C(11)-H(11B) 

H(llA)-  C(11)-H(11B) 

C(11)-C(12)-H(12A) 

C(11)-C(12)-H(12B) 

C(13)-C(12)-H(12A) 

C(13)-C(12)-H(12B) 

H(12A)-  C(12)-H(12B) 

C(8)  -  C(14)  -  H(14) 

C(13)-C(14)-H(14) 

C(15)-C(14)-H(14) 

C(14)-C(15)-H(15A) 

C(14)-C(15)-H(15B) 

C(16)-C(15)-H(15A) 


110.5(2) 

108.4(2) 

112.6(2) 

103.8(2) 

116.8(1) 

109.7(2) 

119.6(2) 

119.2(2) 

118.0(2) 

118.7(2) 

118.1(2) 

122.4(2) 

122.6(2) 

118.3(2) 

119.1(2) 

120.2(2) 

119.9(2) 

119.2(2) 

111.7(1) 

110.5(2) 

106.0(1) 

108.0(1) 

109.1(1) 

112.5(2) 

107.1(2) 

106.4(2) 

109.5(2) 

105.4(2) 

109.0(2) 

111.2(2) 

110.6(2) 

109.1(2) 

108.2(1) 

109.1(1) 

109.7(2) 

112.4(2) 

109.1(1) 

104.2(1) 

108.8(1) 

110.2(1) 

106.9(1) 

109.6(1) 

111.5(2) 

107.3(2) 

107.5(1) 

106.6(1) 

104.4(2) 

106.1(2) 

105.9(1) 

106.8(1) 

108.1(2) 

109.5(2) 

110.0(2) 

111.0(2) 

109.3(2) 

104.5(2) 

110.9(2) 

110.2(2) 

106.5(1) 

108.6(1) 

109.8(2) 

108.3(2) 

110.0(2) 

106.3(2) 

108.5(1) 

110.9(1) 

111.1(1) 

111.9(1) 

105.8(2) 

108.6(2) 

101.8(2) 

107.4(2) 

107.5(1) 

106.0(1) 

109.1(1) 

106.9(1) 

112.1(1) 

114.3(1) 

114.1(1) 

111.5(1) 

111.4(2) 

111.2(2) 

41 


C(16)-C(15)-H(15B) 
H(15A)-  C(15)-H(15B) 
0(3)-  C(16)-H(16) 
C(15)-C(16)-H(16) 
C(17)-C(16)-H(16) 
0(2)-  C(17)-H(17) 
C(13)-C(17)-H(17) 
C(16)-C(17)-H(17) 
C(13)-C(18)-H(18A) 
C(13)-C(18)-H(18B) 
C(13)-C(18)-H(18C) 
H(18A)-  C(18)-H(18B) 
H(18A)-  C(18)-H(18C) 
H(18B)-  C(18)-H(18C) 


109.7(2) 

113.7(2) 

105.8(1) 

102.3(2) 

105.5(1) 

108.3(2) 

111.9(2) 

110.3(2) 

109.1(1) 

106.3(1) 

107.9(2) 

111.0(2) 

106.0(2) 

108.3(2) 

105.7(1) 

106.9(1) 

113.5(1) 

110.5(1) 

107.9(1) 

108.6(1) 

109.8(1) 

113.1(1) 

105.6(1) 

110.7(1) 

109.4(2) 

107.4(2) 

110.5(3) 

106.4(2) 
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Table  5.  Selected  torsion  angles  (°) 


bonds 

Molecule  1 

Molecule  2 

A-ring 

C(10)-C(l)-C(2)-C(3) 

-0.5(2) 

-0-5(2) 

C(l)-C(2)-C(3)-C(4) 

2.1(2) 

0.8(2) 

C(2)-C(3)-C(4)-C(5) 

-2.1(2) 

C(3)-C(4)-C(5)-C(10) 

C(4)-C(5)-C(10)-C(l) 

IIEiH 

C(5)-C(10)-C(l)-C(2) 

-0.3(2) 

B-ring 

C(10)-C(5)-C(6)-C(7) 

-10.5(2) 

C(5)-C(6)-C(7)-C(8) 

41.1(2) _ 

C(6)-C(7)-C(8)-C(9) 

Ewm 

-66.0(2) 

C(7)-C(8)-C(9)-C(10) 

C(8)-C(9)-C(10)-C(5) 

C(9)-C(10)-C(5)-C(6) 

4.2(2) 

C-ring 

C(8)-C(9)-C(ll)-C(12) 

C(9)-C(ll)-C(12)-C(13) 

C(ll)-C(12)-C(13)-C(14) 

UUUM 

56.8(2) 

C(12)-C(13)-C(14)-C(8) 

-60.2(2) 

C(13)-C(14)-C(8)-C(9) 

56.9(2) 

mmsm 

-51.4(2) 

HUDDSHHI 

-31.2(2) 

■om 

2.0(2) 

C(15)-C(16)-C(17)-C(13) 

C(16)-C(17)-C(13)-C(14) 

ESSSH 

C(17)-C(13)-C(14)-C(15) 

Mmsm 

47.7(2) 

C(13)-C(17)-C(16)-0(3) 

150.4(2) 

151.8(2) 

C(15)-C(16)-C(17)-0(2) 

157.7(2) 

152.9(2) 

C(13)-C(17)-0(2)-H(2A) 

-4-7(2) _ 

-78.2(2) 

C(15)-C(16)-0(3)-H(3A) 

UJiUM 

50.8(2) 

C(2)-C(3)-0(1)-H(1A) 

-171.3(2) 

C(4)-C(3)-0(1)-H(1A) 

-8.5(2) 

3.6(2) 
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Table  6.  Hydrogen  bonds 


Bonds 

O-H 

fAl 

o  o 
[Al 

H  O 
fA] 

O-H  O 
Anglen 

molecule  1 

0(1)-H(1A)-  •  •  0(3)' 

0.97 

2.646(2) 

1.682(2) 

171.7(1) 

0(2)-H(2A)-  •  •  0(1)^ 

0.97 

2.773(2) 

1.851(1) 

158.6(1) 

0(3)-H(3A>  .  •  0(1’)^ 

0.97 

2.890(1) 

1.930(1) 

170.5(1) 

1.684(1) 

■I 

yTi«nH:riiJM^B8BBg 

0.97 

0(3’)-H(3’A)-  •  •  0(2)^ 

0.97 

2.736(1) 

1.857(1) 

149.9(1) 

Symmetry  codes:  (l)~x+2,  y-1/2,  -z+1;  (2)  -x+1,  y+1/2,  -z+1;  (3)  x,  y,  z-1  ; 

(4)  -x+1,  y+1/2, -z+2;  (5) -x+2,  y-I/2, -z+2;  (6)  -x+ 1,  y-1/2, -z+1 
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Table  7.  Analysis  of  the  properties  of  the  bond  critical  points 
for  molecule  1  and  2 


Bond 

X, 

^2 

h 

Pc  (eA'^) 

V'pc  (eA-’) 

Rii(A) 

d,(A) 

d2(A) 

£ 

A-ring 

C(l)-C(2) 

-14.91 

-13.13 

9.86 

2,14m 

-18.2(1) 

1.3941 

0.7093 

0.6848 

0.14 

-16.56 

■Elkl 

8.59 

WHtBiW 

-22.2(1) 

1.3968 

0.6823 

0.7145 

0.16 

Bldtlil 

BIcXdel 

i:W 

Hom 

-21.2(1) 

1.3944 

0.7586 

0.6357 

■arei 

BBCT 

MiWiM 

-18.3(1) 

1.3993 

0.6995 

0.6997 

oa 

BEKil 

!>*/<■ 

-17.1(1) 

1.4085 

0.7184 

0.6901 

EMM 

-13.90 

BHETai 

■TiTfla 

HiEiua 

B^m 

1.4064 

0.7400 

0.6664 

iiW 

-15.76 

BEEAl 

wnam 

■liiua 

■Eom 

1.0810 

0.7043 

0.3767 

■«i)rll 

EBEn 

BItfckI 

■wsa 

■Miua 

HBHB 

^^1 

wmm 

0.3943 

■tfia 

WSIBi 

mmfM 

lUiB 

■ftMcia 

■EEHH 

■ilia 

-18.43 

-21.9(9) 

■KBSBI 

■iMtHa 

wsn^ 

-16.37 

11.96 

2.15(2) 

1.3785 

0.8625 

0.5160 

liiUflilUtlH 

gHKl 

-41.68 

28.43 

2.06(8) 

i^HB 

0.9705 

0.8183 

0.1523 

■ifiV 

-39.51 

-38.86 

33.10 

2.19(8) 

e^h 

0.9656 

0.7918 

0.1738 

■>iia 

B-ring 

C(5)-C(6) 

-11.17 

RB 

-11-3(1) _ 

1.5150 

0.7758 

0.7392 

0.05 

iggi 

BlilCT 

lliilB 

■Rsm 

BEsam 

■EBgal 

EHH 

■IBM 

DEH 

WilEB 

BMItB 

■^aa 

BHH 

mam 

■IB1« 

-11.24 

9.55 

■B^ua 

EEHEH 

1.5449 

0.7774 

0.7675 

■ataai 

■UUJi 

wacEB 

■tana 

■Hiiua 

EB9EBi 

1.5278 

0.7819 

0.7459 

■iltiB 

BldHrl 

BE^ 

iHIrl 

■Buoa 

e^ieh 

1.0949 

0.6650 

0.4299 

ww 

Blsilsl 

-13.26 

il:m 

■rtEfeta 

1.0919 

0.7474 

■imEGB 

■IBM 

■slUBlIiE&iH 

-17.45 

-16.39 

■Bueia 

E^m 

1.0934 

0.7486 

■iltEEM 

C(7)-H(7B) 

-15.11 

■Hia 

■BMIEia 

b^dbi 

1.0918 

0.7198 

0.3721 

■IBK# 

Biiita 

■feMcia 

-16.7(2) 

EHa 

■tltVEM 

■tftHI 

Boa 

BHEH 

■BHBa 

-16.0(1) 

Bl^l 

liliEB 

-10.88 

RH 

RB 

Hi 

BSRl 

1.5392 

0.7584 

0.06 

C(I1)-C(12) 

-11.08 

itiiiti 

1-64(1) 

-10-4(1) _ 

1.5409 

0.7602 

VSfM 

■BHI 

-10.91 

9.88 

DERn 

EEEmBi 

1.5326 

0.7590 

0.7737 

■iiE^ 

BIW»< 

-10.40 

9.44 

1.61(1) 

-12.2(1) 

1.5511 

0.7811 

0.7701 

EEB 

KIEISiSHI 

BFHB 

-11.67 

9.11 

1-73(1) 

-15.4(1) 

1.5240 

0.7471 

0.7769 

■Mlf 

■ViHliHIMl 

■Kia 

bebi 

IBB 

HBH 

EEim 

1.0940 

0.7254 

WmmM 

kMI 

MUiBifUka 

-16.42 

UilKH 

aifi.itf 

1.0910 

0.7252 

W(kWm 

■IBM 

C(12)-H(12A) 

-15.29 

■E^ 

■iBn 

-13.9(2) 

0.7609 

■IlCKftlgM 

C(12)-H(12B) 

msm 

-14.39 

■EKB 

■HEH 

-16.4(1) 

0.6880 

■ana 

-17.00 

-16.69 

15.47 

-18.2(2) 

1.1014 

0.7906 

0.3108 

ESa 

IsUEIMeH 

-10.55 

-10.14 

9.96 

-107(1) _ 

1.5407 

0.7822 

0.7585 

0,04 

msm 

-11.26 

9.30 

BEERH 

1.5545 

0.7425 

-11.96 

9.60 

hhb 

EEsan 

1.5415 

0.7855 

EEB 

-11.07 

-10.22 

9.64 

■^DB 

-11.7(1) 

1.5388 

0.7186 

0.8202 

C(15)-H(15A) 

■BEB 

■ma 

12.88 

iliUEIfl 

-20.1(1) 

Mt«aa 

BqBall 

0.3614 

■ilfcB 

C(15VH(15B) 

16.09 

■gtfertia 

-14.9(2) 

mtwMnn 

0.3229 

■iBi.ii 

C(16)-H(16) 

-19.64 

BItXtU 

BItMB 

i^afl 

MMSSm 

1.0963 

uowriai 

0.3182 

C(17)-H(17) 

-17.60 

■BiEI 

IIS9 

MMSM 

IHkliUyiHI 

1.1012 

■tlriail 

■ilia 

mmmm 

-11.21 

iBtigg 

■BBH 

IBIKIIf 

1.5369 

■tWKM 

■IBM 

C(18)-H(18A) 

isdiii^s 

Bsa 

16.28 

1.88(5) 

-16.3(2) 

1.0637 

0.7103 

0.3533 

■IBM 

C(18)-H(18B) 

mmtm 

B^ 

13.99 

1.88(5) 

-18.2(2) 

1.0542 

0.6786 

0.3756 

■ilia 

C(18)-H(18C) 

-16.25 

1  -14.78 

■IM3B 

IBMHB 

1.0574 

0.6869 

0.3706 

■iitii 

C(l7)-0(2) 

-14.87 

IBEItW 

BKisg 

■EHI 

IB^SH 

1.4309 

0.8536 

■iIiTiB 

-14.13 

WltftH 

14.41 

1.93(3) 

-12.8(6) 

1.4168 

0.8365 

IBk^ 

-44.89 

24.94 

2.26(6) 

-67.5(7) 

0.9664 

0.8096 

0.1569 

0.06 

lElsEkl 

-42.90 

24.61 

2.36(8) 

-63.8(7) 

0.9684 

0.7946 

0.1738 

0.06 

-14.88 

13.82 

-15.1(7) 

1.4387 

0.8574 

0.5813 

0.06 

MRiB 

■EUI 

-12.85 

14.64 

laniif 

1.4296 

0.8351 

0.5945 

0.10 

■eUa 

-34.31 

42.26 

2.36(7) 

-27.4(5) 

0.9696 

0.7579 

0.2118 

0.03 

1  . 

1  -40.78 

-39.22 

34.73 

2.25(7) 

-45.3(6) 

0.9678 

0.7921 

0.1757 

0.04 

X|,  ^2  and  Xj  are  the  eigenvalues  of  Hessian  matix.  Pc  -electron  density;  V^pc  -  Laplacian  at  the  CP;  di  &  d: 
are  the  distances  from  the  critical  point;  Ry  -  Sum  of  d|  &  d2 ;  e-  bond  ellipticity. 
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Table  8.  Topological  characteristics  of  the  electron  density  at  the  hydrogen 
Bond  critical  point. 


Bond 

pc(eA'2) 

v2pe(eA-') 

Rii(A) 

d,(A) 

d2(A) 

8 

molecule  1 

0(1)-H(1A)-  •  •  0(3)' 

0.23(5) 

5.01(6) 

1.6854 

1.1482 

0.5372 

0.17 

0(2)-H(2A)-  •  •  0(1)2 

0.11(2) 

3.00(1) 

1.9290 

1.2552 

0.6738 

0.82 

0(3)-H(3A)-.-0(1’)2 

0.20(3) 

2.93(1) 

1.9308 

1.2247 

0.7061 

0.17 

molecule  2 

0(l’)-H(rA)- .  •  0(3’)^ 

0.23(5) 

4.96(5) 

1.6861 

1.1194 

0.5668 

0.20 

0(2’)-H(2’A)-  •  •  0(1’)^ 

0.05(2) 

2.24(1) 

2.0792 

1.3508 

0.7284 

1.42 

0(3’)-H(3’A)-  •  •  0(2’)^ 

0.20(2) 

3.33(2) 

1.8746 

1.2160 

0.6586 

0.15 

Symmetry  codes:  (l)-x+2.  y- 1/2, -z+l;  (2)  -x+1,  y+I/2, -z+l;  (3)x,y,z-l  ; 

(4)  -x+l,y+l/2,-z+2;  (5) -x+2,  y-1/2, -z+2;  (6)  -x+1,  y-I/2, -z+l 
Xu  A,2and  X3  are  the  eigenvalues  of  Hessian  matix.  p,.  -electron  density;  V^p,.  -  Laplacian  at  the  CP;  di  &  d2 
are  the  distances  from  the  critical  point;  Ry  -  Sum  of  d,  &  d2 ;  e-  bond  ellipticity. 
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Table  9.  Net  atomic  charges  (e) 


Atom 

charge 

Atom 

charge 

C(l) 

-0.21(3) 

H(l) 

0.16(2) 

C(2) 

-0.06(3) 

H(2) 

0.29(2) 

C(3) 

0.16(3) 

H(4) 

0.31(2) 

C(4) 

-0.06(3) 

H(6A) 

0.20(2) 

C(5) 

-0.00 

H(6B) 

0.12(2) 

C(10) 

-0.04(2) 

H(7A) 

0.20(2) 

C(6) 

-0.36(3) 

H(7B) 

0.24(2) 

C(7) 

-0.48(3) 

H(8) 

0.16(2) 

C(8) 

-0.12(3) 

H(9) 

0.20(2) 

C(9) 

-0.29(3) 

H(llA) 

0.20(2) 

C(ll) 

-0.31(3) 

H(llB) 

0.22(2) 

C(12) 

-0.47(3) 

H(12A) 

0.08(2) 

C(13) 

-0.00 

H(12B) 

0.26(2) 

C(14) 

-0.07(3) 

H(14) 

0.17(2) 

C(15) 

-0.24(3) 

H(15A) 

0.30(2) 

C(16) 

0.05(3) 

H(15B) 

0.21(2) 

C(17) 

-0.08(3) 

H(16) 

0.21(2) 

C(18) 

-0.61(3) 

H(17) 

0.16(2) 

H(18A) 

0.20(2) 

H(18B) 

0.27(2) 

H(18C) 

0.19(2) 

0(1) 

-0.63(3) 

H(1A) 

0.43(3) 

-0.79(3) 

0.35(3) 

0(2) 

-0.86(3) 

H(2A) 

0.48(2) 

-0.87(3) 

0.46(2) 

0(3) 

-0.84(3) 

H(3A) 

0.25(2) 

-0.77(3) 

0.46(2) 

Figure  1 .  ORTEP  view  of  two  molecules  in  the  asymmetric 
unit  showing  75%  probability  thermal  ellipsoids. 


48 


Figure  2a.  Residual  electron  density  (6presid=pexper  -  pmuit)  niap  for  the  molecule  1.  Contour 
interval  is  0.05  eA'^.  The  solid  lines  show  positive  values,  the  negative  contours  are  dashed  lines, 
the  zaro  one  is  dashed. 
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Figure  2b.  Residual  electron  density  (5presid=pexper  -  pmuit)  map  for  the  molecule  2.  Contour 
interval  is  0.05  eA‘^.  The  solid  lines  show  positive  values,  the  negative  contours  are  dashed  lines, 


the  zaro  one  is  dashed. 


Figure  3.  The  multipole  deformation  eleetron  density  (6pmuit=Pmuit  -  Psph)  map 
ring  plane.  Contours  are  the  same  as  in  Fig.  2. 


Figure  4a.  Electrostatic  potential  of  molecule  (1)  in  the  plane  of  0(1),  C(5)  &  C(16)  atoms. 
Contour  interval  is  0.05  eA'*.  Positive  contours  are  solid  lines,  negative  contours  are  dotted  lines, 
and  zero  contour  is  dashed. 
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Figure  4b.  Electrostatic  potential  of  molecule  (2)  in  the  plane  of  €(!'),  C(5')  &  C(16').  Contour 
are  the  same  as  in  Fig4a. 


Molecule  (2) 


Figure  5.  The  electrostatic  potential  in  the  hydrogen  bonding  environment  of  oxygen 
atoms  of  molecule  (1)  &  (2).  Contour  interval  O.OSe/A;  Positive  solid  line,  negative 
dashed  line,  zero  contour  as  dotted  lines. 


Figure  6.  Isosurface  representation  of  electrostatic  Potential  of  molecule  1  (right) 
and  molecule  2  (left).  Surfaces  are  drawn  at  constant  positive  and  negative 
potentials  of  0.5  e/A  (red)  and  -0.1  e/A  (green).  Gray  dots  are  zero  surface. 
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Appendix  C. 


Electron  density  study  of  estrone 
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Electron  Density  Study  of  Estrone 

by 

Nan  Wu,  Elizabeth  A.  Zhurova  and  A.  Alan  Pinkerton 
Department  of  Chemistry,  University  of  Toledo,  Toledo,  OH  43606,  USA 


Abstract 

Many  breast  cancers  are  known  to  be  hormone  dependent,  i.e.  the 
activation  or  inhibition  of  tumor  growth  can  be  influenced  by  estrogens.  Minor  chemical 
changes  in  the  estrogens  structures  are  known  to  cause  fundamental  changes  in  the 
biological  reactions.  The  reason  for  this  is  unknown,  however,  it  has  been  suggested 
(Vander  Kuur,  J.A.;  Hafner,  M.S.;  Christman,  J.K.;  Brooks,  S.C.  Biochemistry,  1993,  32, 
7016-7021)  that  the  behavior  of  different  estrogen  analogs  can  be  related  to  the  physical 
properties  such  as  the  distribution  of  electron  density  and  the  electrostatic  potential.  In 
this  work  we  present  the  first  electron  density  study  of  an  estrone  molecule  based  on  low- 
temperature  X-ray  diffraction  experimental  data.  The  parameterization  of  the  electron 
density  in  terms  of  a  multipole  expansion  allowed  the  direct  space  calculation  of  the 
estrone  electrostatic  potential.  It  exhibits  negative  regions  around  both  oxygen  atoms, 
more  narrow  for  the  O  atom  bonded  to  hydrogen.  There  are  also  significant  negative 
regions  above  and  below  the  aromatic  ring.  These  areas  represent  the  sites  of  the  possible 
electrophilic  attack  or  the  sites  which  match  the  electropositive  regions  of  the  estrogen 
receptor  during  the  activated  complex  formation.  Chemical  bonding  in  the  estrone 
molecule  as  well  as  the  hydrogen  bonding  is  discussed  in  terms  of  Bader's  topological 
theory. 


Introduction 

Estrogens  and  their  derivatives  have  been  extensively  investigated  since  the 
agonistic  effect  of  estradiol  on  most  breast  cancer  tumors  was  determined'.  Estrogens  and 
related  compounds  bind  as  ligands  to  the  estrogen  receptor  (ER)  forming  an  activated 
complex.  A  series  of  events  are  initiated  by  the  ER  complex  resulting  in  the  activation  or 
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repression  of  selective  genes  and  subsequent  induction  or  suppression  of  production  of 
characteristic  proteins.  In  order  to  better  understand  the  role  of  the  ligand  in  ER  mediated 
processes,  extensive  research  into  the  influence  of  the  different  estrogen  derivatives  on 
the  regulation  of  hormone  responsive  genes  has  been  performed  ’ .  Even  though  most 
changes  in  the  estrogen  ligand  affect  its  affinity  for  the  receptor,  this  does  not  necessarily 
correlate  with  the  stimulation  ability  of  the  transcription  of  estrogen  responsive  genes. 
The  mechanism  by  which  these  structurally  altered  ligands  regulate  the  differential 
induction  of  gene  expression  is  currently  unknown.  However,  it  has  been  suggested  that 
the  differences  in  electrostatic  potential  generated  by  the  lone  pair  electrons  of  the 
phenolic  oxygen  and  n-  electrons  of  the  aromatic  ring  are  responsible  for  the  variation  in 
the  regulation  of  hormone  dependent  genes^.  It  is  obvious,  that  the  binding  ability  of  the 
ligand  to  the  receptor  is  not  the  only  important  factor  in  the  biological  activity  of  a 
potential  drug,  secondary  interactions  and  several  domains  of  the  receptor  are  involved. 
But  the  initial  driving  force  towards  binding  to  a  receptor  is  due  to  the  match  between 
topography  of  the  electrostatic  potential  (ESP)  of  the  pharmaceutical  molecule  and  that  of 
the  binding  site. 

In  principle,  the  ESP  distributions  can  be  calculated  theoretically.  But  the 
calculations  of  reliable  ESP’s  from  high  level  ab  initio  methods  are  only  possible  for 
relatively  small  molecules  because  extended  basis  sets  are  required  in  order  to  obtain 
meaningful  results.  A  recent  approach^  has  to  be  mentioned,  which  uses  large  basis  sets 
for  small  fragments  which  in  turn  are  used  as  building  blocks  for  larger  molecules.  In 
spite  of  the  computational  advantage  of  this  method,  it  can  be  only  a  rough 
approximation  to  the  real  object. 

A  better  approach  is  to  map  the  ESP  experimentally  by  interpretation  of  single 
crystal  X-ray  diffraction  data.  Until  recently,  this  type  of  measurement  was  very  time 
consuming,  even  for  smaller  molecules.  The  reeently  developed  new  generation  of  CCD 
Platform  Diffraetometers  with  an  appropriate  cooling  system  allow  us  to  measure  a 
precise  complete  set  of  data  in  2-3  weeks,  even  for  such  large  molecules,  as  estrogens. 
The  developed  multipole  method*"^  allows  the  parameterization  of  the  electron  density  in 
a  crystal  and  a  subsequent  calculation  of  the  ESP  distribution’®. 
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Another  novel  approach  that  can  be  used  in  the  analysis  of  the  electron  density 
(ED)  and  related  properties  is  Bader’s  topological  theory".  This  theory  is  based  on 
quantum  mechanics  and  describes  a  crystal  in  terms  of  the  ED  p(r),  its  gradient  vector 
field  Vp(r),  ED  curvature,  critical  point  (CP)  positions  and  their  characteristics'^.  The 
CP,  the  point  at  which  Vp(rc)=0,  is  determined  by  it’s  rank,  the  number  of  non-zero 
eigenvalues  of  the  diagonalized  Hessian  (or  curvature)  matrix  X,,  and  signature,  the  sum 
of  the  algebraic  signs  of  X,.  The  ED  exhibits  four  kinds  of  non-degenerate  CP’s  of  rank  3: 
maxima  (3,-3),  minima  (3,+3)  and  two  types  of  saddle  points,  (3,+l)  and  (3,-1).  They 
correspond  to  nuclear  positions,  cages,  rings  and  bonds  in  the  crystal  or  molecule.  In  the 
Vp(r)  field,  there  are  pairs  of  gradient  lines  that  originate  at  a  (3,-1)  CP  and  terminate  at 
two  neighboring  nuclei.  They  are  determined  by  the  eigenvector  corresponding  to  the 
only  positive  eigenvalue  of  the  Hessian  A,3at  this  point  and  they  form  the  atomic 
interaction  lines  along  which  the  ED  decreases  for  any  lateral  displacements.  At 
equilibrium,  this  line  is  named  as  a  bond  path  and  the  associated  (3,-1)  point  is  termed  a 
bond  CP.  The  negative  eigenvalues  and  X2  at  the  bond  CP  correspond  to  the  directions 
normal  to  the  bond  path  and  measure  the  degree  of  ED  contraction  towards  this  point; 
X3  >0  measures  the  degree  of  ED  contraction  towards  each  of  the  neighboring  nuclei.  The 
spatial  distribution  of  the  Laplacian  of  the  ED,  V^p(r)  characterizes  the  concentration  and 
depletion  of  electrons  in  each  point  r  of  the  system.  The  sign  of  the  Laplacian  of  the  ED 
at  the  bond  (3,-1)  CP  depends  on  the  relation  between  the  principal  curvatures  of  the  ED 
at  the  CP  and,  therefore  reflects  the  character  of  the  atomic  interactions.  If  the  electrons 
are  locally  concentrated  in  the  bond  CP  (V^p(rc)<0),  then  ED  is  shared  by  both  nuclei 
(shared  interactions).  This  is  typical  for  the  covalent  type  of  bonds.  If  electrons  are 
concentrated  in  each  of  the  atomic  basins  separately  (  V^p(rc)>0  ),  the  interaction  type  is 
closed-shell,  and  it  is  typical  of  ionic  and  hydrogen  bonds,  as  well  as  van  der  Waals 
interactions'^.  Each  bond  path  is  homeomorphically  mirrored  by  a  virial  path,  a  line  of 
maximally  negative  potential  energy  density  linking  the  same  nuclei.  The  presence  of  a 
bond  path  and  its  associated  virial  path  provide  a  'universal  indicator  of  bonding  between 
the  atoms  so  linked'. 
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At  this  early  state  of  experimental  ED  investigations  on  estrogens,  it  is  very 
difficult  to  formulate  a  precise  hypothesis.  However,  we  believe,  we  can  derive 
parameters  from  experimental  data  which  will  identify  physical  properties  of  regions  of  a 
ligand  that  are  related  to  certain  biological  functions.  This  is  the  information  needed  in 
order  to  intelligently  design  therapeutically  improved  drugs.  The  ED  investigation  of 
estrone,  one  of  the  natural  estrogens,  is  a  very  first  step  towards  the  solution  of  this  task. 

Experimental 

A  colorless  rectangular  crystal  (Table  1)  was  re-crystallized  from  an  ethanol: 
ethyl  acetate  solvent  system.  It  was  mounted  on  a  0.1  mm  capillary  with  epoxy  resin  and 
slowly  cooled  dovm  to  92K  with  an  Oxford  Cryostream.  The  X-ray  diffraction 
experiment  was  performed  with  a  Bruker  platform  diffractometer  with  a  2K  CCD 
detector  using  Ag-Ka  radiation.  0.2°  ©  -scans  were  performed  at  a  detector  distance  of  6 
cm.  Two  detector  settings  of  26=  -20  and  -65  and  4  different  cp  settings  were  used.  The 
exposure  time  of  88  sec  for  the  low-angle  data  and  180  sec  for  the  high-angle  data 
allowed  the  use  of  the  full  dynamic  range  of  the  detector.  The  first  100  frames  were 
repeated  at  the  end  of  the  measurement  to  check  the  X-ray  intensity  decay.  Other 
experimental  details  are  in  Table  1 . 

The  integration  was  performed  with  the  program  SAINT An  empirically  chosen 
integration  box  size  of  1 .8x  1 .8x0.7°  for  low-angle  data  and  1.5x1 .5x0.7°  for  high-angle 
data  and  the  profile  fitting  procedure  based  on  strong  (I  >  IOct)  reflections  allowed  us  to 
obtain  the  best  internal  consistency  (Rint=0.0478)  and  preliminary  refinement  R  factors. 
After  the  integration  process,  the  unit  cell  parameters  were  refined  with  all  reflections 
with  sin0/X,  <  0.84  A”'  since  all  of  the  high-angle  reflections  were  relatively  weak.  The 
integration  process  was  repeated  twice  until  convergence  in  cell  parameters  was  obtained. 
An  empirical  absorption  correction  with  the  program  SADABS'^  was  performed.  The 
data  were  averaged  with  the  program  SORTAV 

Refinements  and  results 
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The  crystal  structure  was  solved  by  direct  methods  and  a  preliminary  least- 
squares  refinement  was  carried  out  with  the  SHELXTL  program  suite'*.  The  positions  of 
hydrogen  atoms  were  determined  from  the  difference  Fourier  peaks.  Anisotropic  thermal 
motion  was  considered  for  all  non-hydrogen  atoms,  and  the  latter  were  refined  as 
isotropic.  The  estrone  molecule  with  50%  probability  thermal  ellipsoids  for  non¬ 
hydrogen  atoms  is  shown  in  Fig. la.  The  atomic  coordinates,  atomic  thermal  parameters 
and  the  scale  factor  were  refined  at  this  stage  with  all  data.  Then,  the  so-called  Hansen- 
Coppens  multipole  refinement*  using  the  XD  program  package’  was  performed  with  the 
low-angle  data  (sin0/A,<  1.0  A'').  In  this  model  the  ED  is  approximated  as  a  sum  of  the 
core  and  valence  pseudo-atomic  electron  densities  in  the  form 

P atomic  (••)  =  Pcorc  C'')  +  P,alence  ('^''■)  +  E  (^'''■)  E  ^ 

/=1  m=-l 

Here  Py  and  P/,„  are  monopole  and  higher  multipole  populations;  Ri  are  normalized 
Slater-type  radial  functions;  yim  are  real  spherical  harmonic  angular  functions;  k'  and  k" 
are  the  valence  shell  contraction\expansion  parameters.  Hartree-Fock  densities'’  are  used 
for  the  spherically  averaged  core  and  valence  shells.  The  drawing  of  the  local  multipole 
coordinate  systems  for  each  atom  has  been  deposited  (Fig.  Dl).  For  the  hydrogen  atoms 
the  x-axis  was  chosen  to  be  directed  towards  neighboring  carbon  or  oxygen  atom.  Only 
dipole  population  along  this  direction  was  refined  for  hydrogen  atoms,  other  populations 
were  set  to  zero.  Since  the  molecule  is  rather  large  and  the  crystallographic  space  group 
(P2\2]2])  is  non-centrosymmetric,  the  X-ray  phase  problem  arises^’.  Therefore,  chemical 
and  symmetry  constraints  were  applied  at  the  initial  stages  of  the  multipole  refinement  as 
follows.  Hydrogen  atoms  attaehed  to  the  same  carbon  were  considered  to  be  chemically 
the  same,  i.e.  having  the  same  monopole  and  dipole  populations  and  kappa  sets,  but 
different  thermal  motion.  H(l),  H(2)  and  H(4)  were  also  taken  the  same.  Carbons  C(l)  - 
C(4);  C(5)  &  C(10);  C(6),  C(7),  C(ll)-C(13),  C(15),C(16),  C(18);  C(8),C(9)  &  C(14) 
were  considered  to  be  the  same.  Atoms  C(l)-C(5),  C(10)  &  C(17)  were  considered  to 
have  the  sp^  hybridization  state,  therefore  only  multipoles  directed  same  as  the  <T-bonds 
and  multipoles  perpendicular  to  the  aromatic  ring  were  refined.  An  additional  mirror 
symmetry  constraint  for  the  aromatic  ring  was  imposed.  A  totally  of  8  k'  and  k"  sets 
were  used  (see  Table  2).  k'  and  k”  parameters  were  fixed  to  1.2  and  1.0  correspondingly 
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for  all  hydrogen  atoms.  Reasonable  values  for  the  initial  k'  and  k"  parameters,  close  to 
the  final  ones  were  picked  up  to  get  better  X-ray  phases  in  the  initial  refinement.  The 
multipole  parameters  were  refined  up  to  the  hexadecapole  level  (1=4)  for  all  non¬ 
hydrogen  atoms.  The  molecular  electroneutrality  constraint  was  also  applied.  For 
hydrogen  atoms,  all  the  parameters  including  positional  and  displacement  ones  were 
refined  with  sin  0/X,  <  0.5  A'‘.  The  0-H  and  C-H  bond  lengths  were  fixed  to  the  tabulated 
values^’  according  to  their  hybridization  state. 

The  model,  described  above  was  refined  with  all  reflections  together  with  atomic 
coordinates,  thermal  parameters  and  scale  factor.  Then,  all  the  constraints,  except  those 
for  hydrogen  atoms  attached  to  sp^  carbons  were  removed  gradually  and  the  model  was 
refined  again.  The  highly  correlated  parameters  were  refined  in  the  separate  groups,  so 
that  the  highest  correlation  coefficient  was  0.54.  In  the  final  cycles,  parameters  with 
values  less  than  one  standard  deviation  were  fixed  to  zero.  The  rigid-bond  test  showed 
that  the  differences  of  mean-square  displacement  amplitudes  along  the  interatomic 
vectors  were  less  than  1.3x10'^  A^.  The  resulting  multipole  and  k  parameters  have  been 
deposited  (tables  D1  and  D2).  The  atomic  coordinates,  thermal  parameters  and 
interatomic  distances  are  deposited  as  a  cif  file. 

Then,  in  order  to  get  the  correct  values  of  atomic  charges^^,  all  the  multipole 
parameters  were  set  to  zero,  and  the  so-called  k-  refinement  was  performed.  Only  Py  and 
k’  parameters  together  with  the  scale  factor,  positional  and  thermal  atomic  parameters 
were  refined.  The  resulting  atomic  charges  and  /c’  values  are  in  the  Table  2. 

The  residual  ED  (6presid=  Pexper  -  Pmuit)  and  model  multipole  deformation  ED 
(Spmuit=  Pmuit  -  Psph)  maps  Were  calculated  with  the  low-angle  data  (  sin  0/X,  <  0.9  A"'). 
They  are  presented  in  Figs.  2  and  3. 

The  ESP  is  a  crystal  (or  molecular)  property  directly  related  to  the  ED: 


nr) =2! 


R,-r 


rp(r’) 


</r'. 


(2) 


where  is  the  charge  of  nucleus  A  located  at  R/(,  pit')  is  the  total  ED.  The  partition  of 
the  ED  into  pseudoatomic  fragments  in  the  form  of  a  multipole  expansion  (1)  allows  the 
direct  space  calculation  of  the  ESP  distribution'®,  recently  becoming  popular  in  such 
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crystallographic  studies^"*’^.  The  ESP  for  the  isolated  molecule  taken  from  the  crystal  is 
shown  in  Fig.4. 

The  topological  analysis  of  the  estrone  ED  was  performed.  The  Laplacian  map  is 
reported  in  Fig.5.  The  bond  CP's  in  the  molecule  were  found  as  well  as  those  representing 
the  hydrogen  bonding  in  the  crystal.  They  are  reported  in  Table  3.  The  dipole  moment  of 
the  molecule  in  the  crystal  after  the  multipole  refinement  (10.2(3)  D)  corresponding  to 
the  geometrical  center  of  the  molecule  was  also  calculated. 

Discussion 

The  estrone  crystal  unit  cell  parameters  were  measured  by  Vanden  Bossche^^  and 
structure  was  solved  first  by  Busetta  et  al.^^  in  1973.  In  the  solid  state,  the  estrone 
molecules  are  linked  into  zig-zag  chains  by  hydrogen  bonds  (Fig.  lb).  Though  their 
refinement  was  performed  up  to  R=0.039  (estrone  II  in  their  nomenclature),  all  atomic 
thermal  parameters  P22  were  too  small  compared  with  Pn  and  P33,  which  reflects, 
probably  the  lack  of  the  experimental  data.  Our  unit  cell  parameters  (Table  1)  only 
approximately  match  those  previously  published^^'^^  with  the  largest  difference  being 
0.09  A.  This  is  a  first  electron  density  investigation  of  a  molecule  belonging  to  the  class 
of  estrogens,  to  the  best  of  our  knowledge.  The  molecule  is  rather  large  and  the  space 
group  P2i2i2i  is  non-centrosymmetric.  This  created  some  problems  in  the  least-squares 
refinement  due  to  the  high  correlations  between  the  refined  parameters  and  the  possible 
wrong  reflection  phases.  The  applied  chemical  and  symmetry  constraints  on  the  atomic 
multipole  parameters  (see  above)  allowed  us  to  solve  this  problem^®  and  get  the 
reasonable  results  reported  herein. 

The  atomic  charges  calculated  as  a  difference  between  the  actual  monopole 
population  number  and  the  number  of  outer  shell  electrons  in  the  neutral  atoms  are 
reported  in  Table  2.  They  all  look  reasonable,  being  negative  for  the  oxygens,  positive  for 
the  carbons  next  to  the  oxygens,  negative  for  the  carbons  in  the  aromatic  ring  and  positive 
for  all  hydrogens.  At  the  same  time,  these  numbers  are  far  from  the  formal  charges  one 
could  expect  (for  example,  +1  for  hydrogen  atom),  which,  generally  is  not  a  surprise  for  a 
covalent  crystal.  These  charges  are  consistent  with  those  for  other  organic  molecules 
obtained  by  multipole  analysis  '  . 
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The  residual  electron  density  (the  difference  between  the  experimental  ED  and  the 
calculated  aspherical  multipole  ED)  map  in  the  aromatic  ring  plane  is  presented  in  Fig.2. 
This  map,  as  well  as  the  maps  in  the  other  planes  (not  presented  in  this  paper),  does  not 
contain  any  feature  greater  than  ±  0.1  eA  which  can  be  considered  as  a  random  noise. 
This  is  one  of  the  indications  of  the  quality  of  the  experimental  data  and  least-squares 
fitting  of  the  model.  Fig.3.  contains  the  model  deformation  electron  density  (DED)  maps, 
which  is  a  difference  between  the  calculated  multipole  (aspherical)  ED  and  the 
calculated  ED  of  neutral  spherical  atoms  at  the  same  positions  in  the  crystal.  Each 
expected  chemical  bond  (in  the  plane  of  the  map)  is  represented  by  a  well-shaped  DED 
peak,  though  the  both  lone  pairs  of  01  (Fig.3a)  and  one  of  the  02  atoms  (Fig.3b)  are 
shifted  from  the  plane.  A  slightly  negative  region  was  found  inside  the  aromatic  ring.  The 
polarization  of  the  DED  on  the  HlA-01  bond  towards  01  atom  is  evident  when 
comparing  to  the  DED  on  the  HI -Cl,  H2-C2  and  H4-C4  bonds.  This  also  could  be 
expected^’  since  the  HI  A  atom  is  involved  in  the  hydrogen  bonding  with  the  02  atom  of 
the  neighboring  molecule  (Fig. lb). 

The  electrostatic  potential  isosurfaces  of  the  isolated  estrone  molecule  taken  from 
the  crystal  are  shown  in  Fig.4.  The  positive  isosurface  of  0.5  eA  *  repeats  the  shape  of  the 
estrone  molecule  (green  color).  The  red  areas  represent  the  negative  isosurface  of  -0.2 
eA"'  (Fig.4a)  and  of -0.1  eA''  (Fig.4b).  The  negative  regions  are  located  around  both 
oxygen  atoms,  smaller  for  the  01  atom  bonded  to  the  hydrogen.  The  02  ESP  negative 
area  is  enveloping  the  H12A  atom.  There  are  also  considerable  negative  regions  above 
and  below  the  aromatic  ring,  where  the  n-  electrons  are  expected  to  be.  Tiny  negative 
ESP  regions  were  found  inside  the  aromatic,  C(5)-C(7)-C(10)  and  C(8)-C(12)-C(14) 
rings.  Their  existence  could  be  caused  by  the  chosen  hydrogen  k-  set^"  (k-1.2  and 
k"=1.0),  not  refined  in  this  work.  All  these  negative  ESP  areas  represent  the  sites  of 
possible  electrophilic  attack  or  the  sites  which  match  electropositive  regions  of  the 
estrogen  receptor  during  the  activated  complex  formation. 

The  topological  analysis"  gives  us  a  more  complete  view  of  the  chemical  bonding 
in  the  crystals.  Fig.5  contains  the  Laplacian  map  in  the  aromatic  ring  plane.  It  was 
calculated  for  the  cluster  of  6  neighboring  molecules  to  include  the  hydrogen  bonding. 
The  Laplacian  is  an  important  function  of  ED,  defined  as^'' 
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V^p(r)=5^p(r)/5x^  +  5^p(r)/5y^  +  p(r)/5z^.  (3) 

It  is  invariant  under  a  rotation  of  the  coordinate  system,  and  it  is  equal  to  the  trace  of  the 
Hessian  matrix  H(r)  with  elements  3^p/5xi5xj.  The  Laplacian  is  also  reference-state 
independent,  it  characterizes  the  actual  concentration  (negative  values)  and  depletion 
(positive  values)  of  electrons  in  the  crystal.  In  Fig.5  the  contours  with  the  same  negative 
Laplacian  values  are  shown  as  the  solid  contours  to  make  it  analogous  to  the  DED 
distribution  maps.  For  all  the  carbons  the  concentration  of  electrons  is  mainly  around  the 
atoms  and  along  the  chemical  bonds,  though  for  the  hydrogens  it  is  almost  spherical.  The 
local  ED  depletions  along  the  chemical  bonds  and  the  contractions  in  the  perpendicular 
directions  indicates  the  positions  of  the  bond  critical  points.  For  the  oxygens  this 
concentrations  is  also  in  the  non-bonding  directions,  which  is  a  quite  common  ‘  and 
associated  with  lone  pairs.  Some  ED  depletion  areas  are  found  inside  the  aromatic  ring 
and  around  it. 

The  positions  of  the  critical  points  of  the  ED  and  their  characteristics  contain  a  lot 
of  information  about  the  character  of  the  chemical  bonding  in  the  crystal.  Table  3 
presents  the  bond  (3,-1)  CP's  found  in  the  estrone  molecule  and  the  one  belonging  to  the 
hydrogen  bond.  The  Table  contains  the  values  of  the  ED  p  and  Laplacian  V^p  at  each  CP, 
eigenvalues  A.|,  ^2,  ^3,  interatomic  distances  Ry  and  the  distances  from  the  CP  to  the 
bonding  atoms  di  and  d2.  As  it  might  be  expected,  the  ED  values  at  the  CP  on  the  double 
bond  0(2)-C(17)  and  on  the  aromatic  ring  bonds  are  generally  higher  than  those  for  the 
single  bonds.  The  ED  value  at  the  0(1)-H(1A)  of  2.07  eA’^  is  also  high  compared  to  the 
C-H  bonds.  The  polarization  of  C-0,  0-H  and  C-H  bonds  is  also  evident,  since  the  CP  is 
located  much  closer  to  the  more  electropositive  atom,  i.e.  electrons  are  shifted  towards 
the  more  electronegative  atom.  The  C-C  bonds  are  less  polarized.  A  quite  short  hydrogen 
bond  (1.852  A)  contains  values  of  the  p  and  V^p^^  considered  typical  for  strong  bonds. 
Therefore,  the  hydrogen  bonding  in  the  estrone  crystal  can  not  be  considered  as  a  weak 
one,  as  it  was  supposed  by  Busetta  et  al. 


Conclusion 
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We  have  performed  a  low-temperature  X-ray  diffraction  investigation  of  the 
electron  density  and  related  properties  of  the  estrone  crystal.  The  molecular  electrostatic 
potential  exhibits  the  negative  regions  around  both  oxygen  atoms,  more  narrow  for  the  O 
atom  bonded  to  hydrogen.  There  are  also  significant  negative  regions  above  and  below 
the  aromatic  ring.  The  quantitative  comparison  of  these  regions  as  well  as  the  atomic 
charges  and  the  dipole  moments  between  different  estrogens  with  different  biological 
activities  would  be  necessary  information  for  the  intelligent  design  of  new  drugs.  This 
will  be  the  task  of  our  future  work. 
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Table  1.  Experimental 


Compound  name: 

Chemical  formula: 

Chemical  formula  weight: 
Crystal  density: 

Crystal  size,  mm 
p,  mm‘' 

Unit  cell  parameters: 

a  (A) 
b(A) 
c(A) 

V(A) 

z 

F(OOO) 

Space  group 
Temperature: 

Radiation(A): 
hmin/hmax 
kmin/kmax 
Imin/lmax 
0min/6  max 

Intensity  decay 

Total  No.  of  reflections 

No.  of  independent  reflections 

E-int 

No.  of  reflections  used 
Criterion 

Refinement  on 
Total  No.  of  parameters 


estrone 

C18H22O2 

270.36 

1.263 

0.2x0.25x0.30 

0.05 

7.7618(1) 

9.9536(1) 

18.4055(1) 

1421.97(2) 

4 

584 

P2i2i2i 

92.0(2) 

Ag  K„  (X=0.56086  A) 
-20/20 
-26/25 
-48/48 
1.75/48.26 
0% 

119720 

25945 

0.0478 

8384 

I>2ct  for  sin0/A,<  0.7  A"' 
I>6ct  for  sin0/Z>  0.7  A'' 
F 

810 
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2 

Weighting  scheme  w=l/a 

Multipole  refinement: 

R  0.0224 

Rw  0.0149 

S  0.62 

K-  refinement: 

R  0.0445 

Rw  0.061 1 

S  1.54 
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Table  3.  The  properties  of  bond  critical  points  in  estrone. 


Bond 

p,  eA'^ 

V^p,  eA'^ 

X\ 

X.2 

^3 

Ry,  A 

di,  A 

d2,  A 

0(1)-C(3) 

2.29(1) 

-27.09(5) 

-19.40 

-17.54 

9.85 

1.367 

0.837 

0.530 

0(2)-C(17) 

2.90(1) 

-7.59(7) 

-27.48 

-23.91 

43.80 

1.228 

0.821 

0.407 

C(l)-C(2) 

2.319(8) 

-24.06(3) 

-15.03 

-14.71 

5.68 

1.395 

0.702 

0.693 

C(l)-C(10) 

2.32(1) 

-25.74(3) 

-16.63 

-15.23 

6.12 

1.404 

0.706 

0.698 

C(2)-C(3) 

2.32(1) 

-25.93(3) 

-17.40 

-14.74 

6.21 

1.399 

0.738 

0.661 

C(3)-C(4) 

2.34(1) 

-24.62(4) 

-17.09 

-14.03 

6.50 

1.397 

0.644 

0.753 

C(4)-C(5) 

2.26(1) 

-23.07(1) 

-15.00 

-13.77 

5.70 

1.404 

0.628 

0.776 

C(5)-C(10) 

2.17(1) 

-21.01(4) 

-14.94 

-13.32 

7.25 

1.412 

0.673 

0.739 

C(5)-C(6) 

1.84(1) 

-14.87(3) 

-12.23 

-11.34 

8.70 

1.518 

0.772 

0.746 

C(6)-C(7) 

1.695(8) 

-11.08(2) 

-10.61 

-9.75 

9.28 

1.531 

0.786 

0.745 

C(7)-C(8) 

1.704(8) 

-11.90(2) 

-10.58 

-10.12 

8.80 

1.528 

0.801 

0.727 

C(8)-C(9) 

1.688(9) 

-11.64(3) 

-10.85 

-9.68 

8.89 

1.547 

0.743 

0.804 

C(8)-C(14) 

1.852(9) 

-16.38(3) 

-12.84 

-11.96 

8.43 

1.526 

0.780 

0.746 

C(9)-C(10) 

1.76(1) 

-15.18(3) 

-11.85 

-11.64 

8.31 

1.528 

0.702 

0.825 

C(9)-C(ll) 

1.660(9) 

-11.24(3) 

-10.65 

-9.73 

9.14 

1.541 

0.753 

0.788 

C(ll)-C(12) 

1.672(9) 

-11.25(3) 

-10.50 

-10.05 

9.30 

1.546 

0.758 

0.788 

C(12)-C(13) 

1.780(9) 

-14.79(3) 

-12.11 

-11.19 

8.51 

1.528 

0.792 

0.736 

C(13)-C(14) 

1.70(1) 

-13.63(3) 

-11.03 

-10.75 

8.16 

1.545 

0.756 

0.789 

C(13)-C(17) 

1.82(1) 

-14.34(4) 

-12.17 

-11.49 

9.32 

1.518 

0.780 

0.738 

C(13)-C(18) 

1.722(9) 

-14.30(3) 

-11.64 

-11.04 

8.39 

1.546 

0.731 

0.815 

C(14)-C(15) 

1.705(9) 

-12.56(3) 

-11.29 

-10.26 

8.99 

1.540 

0.767 

0.773 

C(15)-C(16) 

1.667(9) 

-10.27(3) 

-10.23 

-9.76 

9.71 

1.549 

0.754 

0.795 

C(16)-C(17) 

1.80(1) 

-13.92(3) 

-12.03 

-11.72 

9.88 

1.528 

0.766 

0.762 

0(1)-H(1A) 

2.07(1) 

-39.73(3) 

-39.34 

-38.20 

37.81 

0.967 

0.806 

0.161 

C(l)-H(l) 

1.844(6) 

-18.39(2) 

-16.35 

-15.76 

13.73 

1.077 

0.729 

0.348 

C(2)-H(2) 

1.800(6) 

-17.96(2) 

-16.62 

-14.76 

13.42 

1.077 

0.730 

0.347 

C(4)-H(4) 

1.736(7) 

-17.83(2) 

-15.99 

-14.75 

12.91 

1.077 

0.733 

0.344 

C(6)-H(6A) 

1.713(7) 

-18.26(2) 

-15.99 

-15.18 

12.91 

1.092 

0.744 

0.348 

C(6)-H(6B) 

1.682(8) 

-16.91(2) 

-15.19 

-14.90 

13.18 

1.092 

0.745 

0.347 

C(7)-H(7A) 

1.726(8) 

-16.49(2) 

-15.72 

-14.77 

13.99 

1.092 

0.737 

0.355 

C(7)-H(7B) 

1.723(7) 

-16.00(2) 

-15.61 

-14.89 

14.50 

1.092 

0.743 

0.349 
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'C(8)-H(8) 

.  C(9)-H(9) 

C(11)-H(11A) 
C(11)-H(11B) 
C(12)-H(12A) 
C(12)-H(12B) 
C(14)-H(14) 
C(15)-H(15A) 
C(15)-H(15B) 
C(16)-H(16A) 
C(16)-H(16B) 
C(18)-H(18A) 
C(18)-H(18B) 
C(18)-H(18C) 
H(1A)...X2_0(2) 


1.686(7) 

1.746(8) 


1.613(7) 

1.779(8) 


■15.26(2) 

-15.42(3) 


■11.24(2) 

-15.47(3) 


1.690(8) 

1.638(8) 

1.648(9) 

1.778(8) 

1.724(8) 

1.725(8) 

1.689(8) 

1.794(8) 

-14.04(3) 

-12.42(2) 

-12.87(2) 

-17.11(3) 

-16.28(3) 

-14.31(2) 

-13.34(2) 

-15.90(2) 

0.204(2) 

3.143(2) 

14.25  15.33 


0.354 

0.373 

0.356 

0.350 

0.377 

0.374 

0.365 

0.360 

0.347 

0.364 

0.377 

0.351 

0.345 

0.348 

1.238 


Figure  captions 


Fig.l.  a-  The  estrone  molecule  with  50%  probability  thermal  ellipsoids  for  non-hydrogen  atoms,  b-  Packing 
diagram  showing  the  hydrogen  bonds. 

Fig.2.  Residual  electron  density  (6presid=  Pexper  -  Pmuit)  niap  in  the  aromatic  ring  plane.  Contour  interval  is  0.05  eA 
The  solid  lines  show  positive  bpresid  values,  the  negative  contours  are  dotted  lines,  the  zero  one  is  dashed. 

Fig.3.  The  multipole  model  deformation  electron  density  (6pmuit=  pmuit  -  Psph)  maps:  a-  in  the  aromatic  ring  plane;  b- 
in  the  0(2),  C(17),  C(16)  plane.  Contours  are  the  same  as  in  Fig.2. 

Fig.4.  Isosurfaces  of  the  electrostatic  potential  of  the  isolated  estrone  molecule  taken  from  the  crystal.  Green  is  the 
surface  of  0.5  eA"',  red  is  the  one  of -0.2  eA"'  (a)  and  of -0.1  eA  '  (b). 

Fig.5.  The  Laplacian  of  the  electron  density  of  estrone  in  the  aromatic  ring  plane.  The  solid  lines  show  the  negative 
Laplacian  values,  i.e.  they  present  the  accumulation  of  electrons  in  the  crystal.  The  contours  are:  -125,  -110,  -95,  - 
80,  -65,  -50,  -45,  -40,  -35,  -30,  -25,  -20,  -15,  -10,  -5  eA'^  are  solid  lines;  3,  6,  9,  12,  15,  18,  21,  24,  27,  30  eA'^  are 
dotted  lines. 
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Figure  la.  The  estrone  molecule  with  50%  probability  thermal  ellipsoids 
for  non-hydrogen  atoms. 
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Figure  2.  Residual  electron  density  (6presid=Pexper  -  Pmuit)  map  in  the  aromatic  ring  plane.  Contour  interval  is  0.05  eA 
The  solid  lines  show  positive  values,  the  negative  contours  are  dashed  lines,  the  zaro  one  is  dashed. 
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Figure  3a.  The  multipole  deformation  electron  density  (5pmuit=Pmuit-  Psph)  map  in  the  aromatic  ring  plane.  Contours 
are  the  same  as  in  Fig.  2. 
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Figure  3b.  The  multipole  deformation  electron  density  (5pmuii=Pniuit-  Psph)  map  in  the  02,  C17,  C16  plane.  Contours 
are  the  same  as  in  Fig.  2. 
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rr 


Figure  5.  The  Laplacian  of  the  electron  density  of  estrone  in  the  aromatic  ring  plane. 

The  solid  lines  show  the  negative  Laplacian  values,  i.e  they  present  the 
accumulation  of  electrons  in  the  ciystal.  The  contours  are:  -125,  -1 10,  -95, 

-80,  -65,  -50,  -45,  -40,  -35,  -30,  -25,  -20,  -15,  -10,  -5,  3,  6,  9,  12,  15,  18,  21,  24, 
27,  30  eA^. 
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The  deposition  material: 


Fig.Dl. 
Table  Dl. 
Table  D2. 
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Figure  Dl. 


Local  coordinate  system  in  estrone 
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Ted>le  Dl.  The  multipole  population  parameters  for  non-hydrogen  atoms  in  estrone 


Table  Dl.  (continued) 


Table  D2.  The  multipole  population  parameters  for  hydrogen  atoms  in  estrone 


H(llA) 

_ 

0.756(6) 

1.2 

1.0 

0.097(8) 

H(9) 

(i)sro 

01 

Zl 

(8)0080 

H(8) 

_ 

0.729(8) 

1.2 

1.0 

0.12(1) 

H(7B) 

0.716(6) 

1.2 

1.0 

0.156(8) 

H(7A) 

0.716(6) 

1.2 

1.0 

0.156(8) 

H(6B) 

0.658(6) 

1.2 

1.0 

0.190(7) 

H(6A) 

0.658(6) 

1.2 

1.0 

0.190(7) 

*s 

W 

X 

0.664(8) 

1.2 

1.0 

0.18(1) 

S 

X 

0.702(8) 

1.2 

1.0 

0.21(1) 

(l)H 

0.729(8) 

1.2 

1.0 

0.20(1) 

H(1A) 

0.574(8) 

1.2 

1.0 

0.10(1) 

Atom 

i 

oT  y  V  O 

C 

o 

o 


o 

2 

H 


H(18C) 

0.762(5) 

1.2 

1.0 

0.140(7) 

H(18B) 

0.762(5) 

1.2 

1.0 

0.140(7) 

H(18A) 

_ 

0.762(5) 

1.2 

1.0 

0.140(7) 

H(16B) 

0.757(6) 

1.2 

1.0 

0.198(8) 

H(16A) 

0.757(6) 

1.2 

1.0 

0.198(8) 

H(15B) 

0.741(6) 

1.2 

1.0 

0.15(1) 

(VSl)H 

0.741(6) 

1,2 

1.0 

0.15(1) 

H(14) 

0.777(8) 

1.2 

1.0 

0.12(1) 

H(12B) 

0.835(6) 

1.2 

1.0 

0.149(9) 

H(12A) 

0.835(6) 

1.2 

1.0 

0.149(9) 

H(llB) 

0.756(6) 

1.2 

1.0 

0.097(8) 

Atom 

+ 

Cu  *y  Q 

VO 

00 


Appendix  D. 


Crystal  structure  of  4-HydroxyestradioI 
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Tcdale  1 .  Crystal  data  and 

Crystal  Data 

Empirical  formula 
Formula  weight 
Temperature 
Crystal  system 
Space  group 
Unit  cell  dimensions 


Volume,  Z 

Density  (calculated) 
Absorption  coefficient 
F(OOO) 

Crystal  size 
Crystal  color 

Data  Collection 

Diffractometer 
Radiation,  Wavelength 
Scan  method 

Crystal-Detector  distance 

Gmax  for  data  collection 
HKL  Limiting  indices 

^int 

Unique  Data 

Refinement 

No.  of  Reflections 
R 


Experimental  details 


C18H24O3 

288.37 
120  K 

Monoclinic 

P2i 

a  =  5.651(1)A 
P  =  107.394  (3)° 
b  =  19.169(3)A 
c  =  6.980(1)A 
721.5(2)A^  ,2 
1.327  Mg/m^ 

0.08  9  mm"^ 

312 

0 . 35x  0.3  X  0.21  mm 
Pale  yellow 

Brukers'  SMART  IK  CCD 
Mo  Ka  ,  1=0.71073  A 

(O-scan,  Scan  width=  0.3° 

4 . 53cm 

25.37° 

“6  ^ —  h  ^ — 6,  ~22  4 —  k  ^ —  22, 

-8  <-  1  <-7 
0.031 
2530 


2340  [Fo  >4a(Fo)] 
0.044 
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wR 

0.10 

S  (  GoodnesS“Of“fit ) 

1.11 

No.  of  parameters  refined 

286 

(Shift/esd)max 

0.001 

'^Pmax  /^Pmin 

0.18/-0.24  eA'^ 
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Figure  1 .  ORTEP  view  of  4-Hydroxyestradiol  showing  50%  probability  thermal  ellipsoids. 
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Tcible  2 .  Fractional  atomic  coordinates  (xlO”^)  and  equivalent 
isotropic  displacement  Parameters  (A^xlO^)  . 


Atom 

X 

y 

z 

U(eq) 

0(1) 

5626(4) 

826(1) 

1077  (3) 

24(1) 

0(2) 

6651  (3) 

276(1) 

4721  (3) 

22(1) 

0(3) 

-4509(3) 

-3830(1) 

2856(3) 

22(1) 

C(l) 

339(5) 

-384  (1) 

362(4) 

20(1) 

C(2) 

1866(5) 

117(1) 

-122(4) 

20(1) 

C(3) 

3992(5) 

330(1) 

1324  (4) 

20(1) 

C(4) 

4551(5) 

43(1) 

3244  (4) 

18(1) 

C(5) 

3066(5) 

-457(1) 

3750(4) 

18(1) 

C(6) 

3879(5) 

-762(1) 

5854  (4) 

19(1) 

C(7) 

1970 (5) 

-1253(1) 

6276(4) 

21(1) 

C(8) 

769(5) 

-1711(1) 

4465  (4  ) 

17(1) 

C(9) 

-757  (5) 

-1241(1) 

2745  (4) 

18(1) 

C(10) 

900(4) 

-680(1) 

2258 (4) 

17(1) 

C(ll) 

-2263(5) 

-1657(1) 

903(4) 

22(1) 

C(12) 

-3810(5) 

-2244  (2) 

1434(4) 

22(1) 

C(13) 

-2211(4) 

-2720(1) 

3060 (4) 

18(1) 

C(14) 

-888  (5) 

-2273(1) 

4900(4) 

18(1) 

C(15) 

227  (5) 

-2807  (1) 

6573(4) 

23(1) 

C(16) 

-1728(5) 

-3397(1) 

6130(4) 

23(1) 

C(17) 

-3598(5) 

-3225(1) 

4064 (4) 

18(1) 

0(18) 

-399(5) 

-3135(1) 

2239(4) 

21(1) 

Ueq  is  defined  as  one  third  of  the  trace  of  the  orthogonalized 
Uij  tesor. 


Table  3 .  Bond  lengths (A) 


0(1)-C{3) 

0(1)-‘H(1A) 

0(2)>C(4) 

0(2)-H(2A) 

0(3)-C(17) 

0(3) -H(3A) 

C(l)-C{10) 

C(l)-C(2) 

C(l)-H(l) 

C(2)-C(3) 

C(2)-H(2) 

C(3)-C(4) 

C(4)-C{5) 

C(5)-C(10) 

C(5)-C(6) 

C(6)-C(7) 

C(6)“H(6B) 

C(6)“H(6A) 

C(7)-C(8) 

C(7)-H(7B) 

C(7)-H{7A) 

C(8)’’C(14) 

C(8)-C(9) 

C(8)-H(8) 

C(9)-C(10) 

C(9)-C(ll) 

C(9)-H(9) 

C(ll)-C(12) 

C(ll) -H(llB) 

C(11)-H(11A) 

C(12)-C(13) 

C{12) “H(12A) 

C(12)-H(12B) 

C(13) -C(18) 

C(13)-C(14) 


1.371 (3) 

0.86(4) 

1.393(3) 

0.82  (4) 

1.436(3) 

0. 90  (4) 

1.387(3) 

1.398  (4) 

0.98  (2) 

1.380  (4) 

0.95(3) 

1.395(3) 

1.387(4) 

1.416(3) 

1.518 (3) 

1.525(4) 

0.98 (3) 

1.02 (3) 

1.523(4) 

0.99(3) 

0.95 (3) 

1.517 (3) 

1.542  (3) 

0.99(3) 

1.530(3) 

1.536(4) 

0.99(3) 

1.536(4) 

1.00  (3) 

0.95(3) 

1.525(4) 

0.94 (3) 

0.95 (3) 

1.536(4) 

1.539(3) 
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C(13)-C(17) 

1.539(3) 

C(14)-C(15) 

1.538(4) 

C(14)-H(14) 

1.00(3) 

C(15)-C(16) 

1.547(4) 

C(15)-H(15A) 

0.98 (3) 

C(15)-H(15B) 

1.06(3) 

C(16)-C(17) 

1.547(3) 

C(16)-H(16B) 

0.95(3) 

C(16)-H(16A) 

1.01 (3) 

C(17)-H(17) 

1.00(3) 

C(18)-H(18B) 

1.01 (3) 

C(18)-H(18A) 

0.98 (3) 

C(18)-H(18C) 

0.97 (3) 

Bond  angles  (*’) 


C{10)-C(l)-C(2) 

122.0(2) 

C(3)-C(2)-C(l) 

119.5(2) 

0(1)-C(3)-C(2) 

125.8(2) 

0(1)-C(3)-C(4) 

115.3(2) 

C(2)-C(3)-C(4) 

118.9(2) 

C(5)-C(4)-0(2) 

118.6(2) 

C(5)-C(4)-C(3) 

122.5(2) 

0{2)-C(4)-C(3) 

118.9(2) 

C(4)-C(5)-C(10) 

118.5(2) 

C(4)-C(5)-C(6) 

118.8(2) 

C(10)-C(5)-C(6) 

122.7(2) 

C(5)-C(6)-C(7) 

113.4  (2) 

C(8)-C(7)-C(6) 

111.2  (2) 

C(14)-C(8)-C(7) 

113.1(2) 

C(14)-C(8)-C(9) 

109.6(2) 

C(7)-C(8)-C(9) 

108.4(2) 

C(10) -C(9)-C(ll) 

113.7(2) 

C(10)-C(9)-C{8) 

110.3(2) 

C(ll)-C(9)-C(8) 

112.9(2) 

C(l)-C(10)-C(5) 

118.6(2) 

C(l)-C(10)-C(9) 

121.3(2) 

C(5)-C(10)-C(9) 

120.1(2) 

C(12) -C(ll)-C(9) 

113.0(2) 

C(13)-C{12)-C(ll) 

111.4  (2) 

C(12)-C(13)-C(18) 

110.2  (2) 

C(12)-C(13)-C(14) 

108.7(2) 

C(18) -C(13)-C(14) 

112.8(2) 

C(12)-C(13)-C(17) 

116.5  (2) 

C(18)-C(13)-C{17) 

109.3(2) 

C(14)-C(13)-C(17) 

99.0  (2) 

C(8)-C(14)-C(15) 

119.2(2) 

C(8)-C(14)-C(13) 

113.8  (2) 

C(15)-C(14)-C(13) 

104.5(2) 

C(14) -C(15)-C(16) 

103.2(2) 

C(15)-C(16)-C(17) 

106.1(2) 

0(3)-C(17)-C(13) 

112.5  (2) 

0(3)-C(17)-C(16) 

113.6(2) 

C(13)-C(17)-C(16) 

105.0  (2) 

C(3)-0(1)-H(1A) 

108 (2) 

C(4)-0(2)-H(2A) 

104 (2) 

C(17) -0(3)-H(3A) 

105(3) 

C(10)-C(l)-H(l) 

119.1(15) 

C(2)-C(l)-H(l) 

118.8 (15) 

C(3)-C(2)-H(2) 

120.4(15) 

C{l)-C(2)-H(2) 

120.2(15) 

C(5)-C(6) -H(6B) 

108.8(18) 

C{7)-C(6)-H(6B) 

110.3(18) 

C(5)-C(6)-H(6A) 

109.3(16) 

C(7)-C(6)-H(6A) 

110.1(17) 

H(6B)-C(6)-H(6A) 

104 (2) 

C{8)-C(7)-H(7B) 

111.9(17) 

C(6)-C(7) -H(7B) 

108.2(17) 

C(8)-C(7)-H(7A) 

109.4 (17) 

C(6)-C(7) -H(7A) 

112.6(17) 

H(7B)-C(7)-H{7A) 

103  (2) 

C(14)-C(8)-H(8) 

112.6(15) 

C(7)-C(8)-H(8) 

107.4  (15) 

C(9)-C(8)-H(8) 

105.3(15) 

C(10)-C(9)-H(9) 

109.5(14) 

C{ll)-C(9)-H(9) 

106.6(14) 

C(8)-C(9)-H(9) 

103.3(14) 

C(12)-C(11)-H(11B) 

114.8(17) 

C(9)-C(11)-H(11B) 

106.7(16) 

C{12)-C(11)-H(11A) 

109.0(16) 

C(9)-C(11)-H(11A) 

110.3(16) 

H{11B)-C(11)-H(11A) 

103(2) 

C(13)-C(12)-H(12A) 

112.0(19) 

C(11)-C(12)-H(12A) 

112.8(18) 

C{13)-C(12)-H(12B) 

107.3(17) 

C(11)-C(12)-H{12B) 

112.3(18) 

H(12A)-C(12)-H(12B) 

101  (2) 

C(8)-C(14)-H(14) 

105.1(15) 

C(15)-C(14)-H(14) 

108.1(15) 

C(13)-C(14)-H(14) 

10.5.3(15) 

C(14)-C(15)-H(15A) 

111.0(17) 

C(16)-C(15)-H(15A) 

113.6(18) 

C{14)-C(15)-H(15B) 

111.0(16) 

C(16)-C(15)-H(15B) 

115.4(16) 

H(15A)-C(15)-H(15B) 

103(2) 

C(15)-C(16)-H(16B) 

110.4(17) 

C(17)-C(16)-H(16B) 

108.7(16) 

C(15)-C(16)-H(16A) 

111.1(16) 

C(17)-C(16)-H(16A) 

109.5(15) 

H(16B)-C(16)-H(16A) 

111(2) 

0(3)-C(17)-H(17) 

106.7(14) 

C(13)-C(17)-H(17) 

109.9(15) 

C(16)-C(17)-H(17) 

109.1(13) 

C(13)-C(18)-H(18B) 

110.5(17) 

C(13) -C(18) -H (18A) 

111.9(16) 

H(18B)-C(18)-H(18A) 

103(2) 

C(13)-C(18)-H(18C) 

112.0(15) 

H(18B)-C(18)-H(18C) 

112(2) 

H(18A)-C(18)-H(18C) 

107  (2) 

Tcible  4.  Aniscbropic  displacemen-t  parameters  (A  x  10^) 
for  non  hydrogen  atoms . 


Atom 

Ull 

U22 

U33 

U23 

U13 

U12 

0(1) 

27(1) 

20(1) 

25(1) 

2(1) 

9(1) 

-7(1) 

0(2) 

23(1) 

19(1) 

23(1) 

-1(1) 

6(1) 

-6(1) 

0(3) 

26(1) 

18(1) 

22(1) 

-2(1) 

9(1) 

-5(1) 

C(l) 

19(1) 

17(1) 

21(1) 

-2(1) 

4  (1) 

0(1) 

C(2) 

25(1) 

15(1) 

18  (1) 

2(1) 

5(1) 

2(1) 

C(3) 

24  (1) 

13(1) 

27  (1) 

-1(1) 

12(1) 

2(1) 

0(4) 

18(1) 

16(1) 

18  (1) 

-4  (1) 

4  (1) 

0(1) 

0(5) 

21(1) 

14  (1) 

20(1) 

-1(1) 

9(1) 

4  (1) 

0(6) 

20(1) 

17(1) 

18(1) 

-1(1) 

3(1) 

0(1) 

0(7) 

26(2) 

22(1) 

16(1) 

2(1) 

7(1) 

-2(1) 

0(8) 

16(1) 

18(1) 

16(1) 

0(1) 

5(1) 

1(1) 

0(9) 

18(1) 

16(1) 

22(1) 

1(1) 

8(1) 

1(1) 

0(10) 

20(1) 

12(1) 

20(1) 

-3(1) 

7(1) 

1(1) 

0(11) 

21(1) 

24  (1) 

15(1) 

6(1) 

-4(1) 

-5(1) 

0(12) 

19(1) 

22(1) 

22(1) 

-1(1) 

2(1) 

-4(1) 

0(13) 

16(1) 

16(1) 

21(1) 

0(1) 

6(1) 

-1(1) 

0(14) 

19(1) 

17(1) 

18  (1) 

0(1) 

6(1) 

1(1) 

0(15) 

28  (2) 

21(1) 

18  (1) 

0(1) 

3(1) 

-5(1) 

0(16) 

26(1) 

20  (2) 

23(1) 

3(1) 

6(1) 

-2(1) 

0(17) 

20(1) 

15(1) 

22  (1) 

-2(1) 

9(1) 

-1(1) 

0(18) 

22(1) 

19(1) 

22  (1) 

-4(1) 

8(1) 

-5(1) 

The  anisotropic 

-2n^  [h^  a*^Uii+ 

displacement  factor 

.  ..  +2hka*b*Ui2] 

exponent  takes 

the  form: 
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Table  5 .  Hydrogen  coordinates  (A  x  10^)  and  isotropic 
displacement  parameters  (A^  x  10^)  . 


Atom 

X 

y 

z 

Ueq 

H(9) 

-1980 (50) 

-1019(13) 

3320(40) 

13(7) 

H(8) 

2120 (50) 

-1905(13) 

3990(40) 

10(6) 

H  (12A) 

-4740(60) 

-2497(16) 

300(50) 

27  (8) 

H(18B) 

-1320  (50) 

-3379(17) 

950  (40) 

30  (8) 

H(16B) 

-2590(50) 

-3401(15) 

7120(40) 

20  (7) 

H(2) 

1450(50) 

304  (14  ) 

-1440 (40) 

13(6) 

H(15A) 

1870(60) 

-2960(16) 

6540  (40) 

26(8) 

H(6B) 

4250(60) 

-379(17) 

6830(50) 

29  (8) 

H(llB) 

-1070(50) 

-1814  (15) 

190(40) 

19(7) 

H  {16A) 

-930(50) 

-3867(17) 

6060 (40) 

22  (7) 

H  (18A) 

750(50) 

-2827(15) 

1820  (40) 

18(7) 

H(llA) 

-3330  (50) 

-1357(15) 

-40(40) 

14  (6) 

H(7B) 

2810  (50) 

-1532(16) 

7480  (50) 

23  (7) 

H(6A) 

5520(60) 

-1013(16) 

6080(40) 

26(8) 

H(7A) 

720(50) 

-1011(15) 

6660  (40) 

18(7) 

H (15B) 

610(50) 

-2572(16) 

8020  (50) 

31(8) 

H(12B) 

-5110(60) 

-2068(15) 

1910(40) 

18(7) 

H(l) 

-1130  (50) 

-546(13) 

-690(40) 

11(6) 

H(18C) 

610(50) 

-3455(15) 

3230(40) 

15(7) 

H(1A) 

5230(60) 

933(17) 

-170  (50) 

35(9) 

H(14) 

-2240(50) 

-2016(14) 

5250(40) 

16(7) 

H(2A) 

7460(60) 

470  (20) 

4090(50) 

39(10) 

H(3A) 

-5470(70) 

-4060 (20) 

3490(60) 

59(12) 

H(17) 

-5070(50) 

-2987  (14) 

4270(30) 

9(6) 
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Table  6 .  Selected  torsion  angles ( ° ) 


C(10)-C(l)-C(2)-C(3) 

■”0.4(4) 

C(1)-C(2)~C(3)-0(1) 

-178.5(2) 

C(l)-C(2)-C(3)-C(4) 

-0.6(4) 

0(1)-C(3)-C(4)-C(5) 

179.0 (2) 

C(2)-C(3)-C(4)-C{5) 

0.9(4) 

0(l)-C{3)-C(4)-0(2) 

0.3(3) 

C(2)-C{3)-C(4)-0(2) 

-177.8(2) 

O(2)-C(4)-C(5)-C(10) 

178.4  (2) 

C(3)-C(4)-C(5)-C(10) 

-0.2(4) 

0(2)-C(4)-C(5)-C(6) 

-3.8  (3) 

C(3)-C(4)-C(5)-C(6) 

177.5  (2) 

C(4)-C(5)-C(6)-C{7) 

174 .5  (2) 

C(10)-C(5)-C{6)-C{7) 

-7.8(3) 

C(5)’C(6)-C(7)-C(8) 

40.0(3) 

C(6)-C{7)-C(8)-C(14) 

172.3(2) 

C(6)-C(7)-C(8)-C(9) 

-65.9  (3) 

C(14)-C(8)-C{9)-C(10) 

-178.9  (2) 

C(7)-C(8) “C(9)-C{10) 

57.2(3) 

C(14)-C(8)“C{9)-C(11) 

-50.4  (3) 

C(7)-C(8)-C(9)-C(ll) 

-174.3(2) 

C(2)“C{1)-C(10)-C(5) 

1.1(4) 

C(2)“C{1)“C(10)-C(9) 

-178.4  (2) 

C(4)“C(5)-C(10)-C(1) 

-0.7(3) 

C(6)-C(5)-C(10)-C(l) 

-178.4  (2) 

C(4)-C(5)“C(10)-“C(9) 

178.7  (2) 

C(6)-C(5)-C(10)-C(9) 

1.1(4) 

C(11)-C(9)-C{10)-C{1) 

25.5(3) 

C(8)-C(9) -C(IO) -C(l) 

153.5  (2) 

C(ll)-C(9)-C(10)-C(5) 

-154 .0  (2) 

C(8)-C(9)-C(10)-C(5) 

-26.0(3) 

C(10)  -C(9)-C(ll)*-C(12) 

177.0  (2) 

C(8)-C(9)-C(ll)-C(12) 

50.3(3) 

C(9)-C(11)“C{12)-’C(13) 

-53.4 (3) 

C(11)-C(12)“C(13)--C(18) 

-68.0(3) 

C(ll)-C(12)-C(13)-C(14) 

56.1(3) 

C(11)-C(12)-C(13)“C(17) 

166.8  (2) 
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C{7)-C(8)-C(14)-C{15) 

C{9)-C(8)-C(14)-C(15) 

C(7)-C(8)-C(14)-C(13) 

C(9)-C(8)-C(14)-C(13) 

C(12)-C(13)-C{14)-C(8) 

C(18)-C(13)-C(14)-C(8) 

C(17)-C(13)-C(14)-C(8) 

C(12)-C(13)-C(14)-C(15) 

C{18)-C(13)-C(14)-C(15) 

C{17)-C(13)-C(14)-C(15) 

C(8)-C(14)-C(15)-C(16) 

C(13)-C(14)-C{15)-C(16) 

C(14)-C(15)-C{16)-C(17) 

C(12)-C(13)-C(17)-0(3) 

C(18)-C(13)-C(17)-0(3) 

C(14)-C(13)-C(17)-0(3) 

C(12)-C(13)-C{17)-C(16) 

C(18)-C(13)-C(17)-C(16) 

C(14) -C(13)-C(17) -C(16) 

C(15)-C(16)-C(17)-0(3) 

C(15)-C(16)-C(17)-C(13) 


-58.7(3) 
-179.8  (2) 
177.3  (  i 
56.2(3) 
-59.5(3) 
63.0(3) 
178.4(2) 
168.9(2) 
-68.6(3) 
46.9(2) 
-163.5(2) 
-35.1(3) 
9.0(3) 
79.3(3) 
-46.4 (3) 
-164.5(2) 
-156.6(2) 
77.7  (2) 
-40.5(2) 
143.4 (2) 
20.0(3) 
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